r 



FILE COPY 



NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 



REPORT No. 773 

ANALYSIS OF HEAT AND COMPRESSIBILITY EFFECTS 
IN INTERNAL FLOW SYSTEMS AND HIGH-SPEED 
TESTS OF A RAM-JET SYSTEM 



By JOHN V. BECKER and DONALD D. BAALS 



NATION/ L 
LANOLE 



*OUESTS 
^ follows: 



FO 



NATIONAL 

1724 r 
WASHINGTON 



ON LOAN FROM THE FILES Of 




ADVISORY COMMITTEE FOR AERONAUTICS 
MEMORIAL ^ZROmrm LABORATORy 
HELD. HAMPTON, VIRGINIA 



ABOVE ADDRESS. 



PUBLICATIONS SHOULD BE AOORE^<;f_o 



ADVISORY COMMITTEE FOR AEffONA. 
N.W., 
25, D.C. 



SI ^EET. 



1943 



For sale by the Superintendent of Documents, XJ. S. GoTernment Printing Ofiice. Washington 25, D. C. 



Price 30 cents 



AERONAUTIC SYMBOLS 
1. FUNDAMENTAL AND DERIVED UNITS 



Symbol 



Length- 
Time—- 
Force- _. 

Power. _ 
Speed- - 



I 
i 

F 



P 
V 



Metric 



Unit 



meter 

second 

weight of 1 kilogram- 



horsepower (metric)- 
fkilometers per hour. 
\meters per second - _ 



Abbrevia- 
tion 



s 

kg 



kph 
mps 



English 



Unit 



foot (or mile) 

second (or hour)__ 
weight of 1 pound 

horsepower 

miles per hour 

feet per second 



Abbrevia- 
tion 



ft (or mi) 
sec (or hr) 

lb 



mpn 
fps 



2. GENERAL SYMBOLS 



W 
9 

m 
I 



S 

G 

b 

c 

A 
V 

L 
D 
D, 
Di 

Op 
O 



Weight =m5r 

Standard acceleration of gravity=9.80665 m/s^ 
or 32.1740 ft/sec^ 

Mass=— 
9 

Moment of inertia=m^r2. (Indicate axis of 

radius of gyration k by proper subscript.) 
Coefficient of viscosity 



p Kinematic viscosity 

p Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-m'^-s^ at 15° C 

and 760 mm; or 0.002378 Ib-ft"^ sec^ 
Specific weight of '^standard'' air, 1.2255 kg/m^ or 

0.07651 Ib/cu ft 



3. AERODYNAMIC SYMBOLS 



Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio, ^ 
True air speed 
Dynamic pressure, ^pV^ 

Lift, absolute coefficient Cl=-^ 

Drag, absolute coefficient Cd=-^ 

Profile drag, absolute coefficient (7do=^ 

Induced drag, absolute coefficient (7d<=^ 

Parasite drag, absolute coefficient (^dp='^ 

Q 

Cross-wind force, absolute coefficient Cc=-^ 



itc Angle of setting of wings (relative to thrust line) 

it Angle of stabiUzer setting (relative to thrust 
line) 

Q Resultant moment 

12 Resultant angular velocity 

B Reynolds number, p-^ where Z is a linear dimen- 
sion (e.g., for an airfoil of 1 .0 ft chord, 100 mph, 
standard pressiu^e at 15° C, the corresponding 
Reynolds number is 935,400; or for an airfoil 
of 1.0 m chord, 100 mps, the corresponding 
Reynolds number is 6,865,000) 

a Angle of attack 

€ Angle of downwash 

ao Angle of attack, infinite aspect ratio 

at Angle of attack, induced 

tta Angle of attack, absolute (measured from zero- 
lift position) 
7 Fhght-path angle 
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ANALYSIS OF HEAT AND COMPRESSIBILITY EFFECTS IN INTERNAL FLOW SYSTEMS AND 

HIGH-SPEED TESTS OF A RAM-JET SYSTEM 

]3v John \'. Beckkr and Donald ]). Baals 



SUMMARY 

An analysis has been made by the NACA of the effects of heat 
arid compressibility in the flow through the internal systetns of 
aircraft. Equations and charts are developed whereby the flow 
characteristics at key stations in a typical internal system may 
be readily obtained. The effects of compressibility (density 
change) were fully accounted for and were found to be important 
at present-day flight conditions, particularly at high altitudes. 
Particular attention was given to the flow changes across radiators 
and air-cooled engines. It is shown that very high Mach 
numbers will be attained vnthin the fins of existing air-cooled^ 
engines if required cooling is to be maintained at altitude. 
The density decrease and the pressure drop across radiators and 
engines are shown to be considerably greater than the values 
computed on the basis of an incompressible flow. The equations 
show that when a Mach number of 1.0 is reached at any station 
of the internal system, the flow becomes choked and any further 
reduction of outlet pressvre will not increase the floio rate. 
The addition of heat had a throttling effect on the flow rate and 
led to choking at lower rates of flow. 

Tests of a 13-inch didnidrr ram-jet type of propulsion system 
using an electrical heaiing decice to add heat to the internal flow 

% - > • «' T. 1 O— C 1 T • I > ' - ' 
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('oiiipr(4ionsiv(» analysis employing compressiblc-llow rela- 
tions, tlioreforo, is roquirod in ordcM* to compute accurately 
the conditions at each station within the duct, the pressure 
drop across the radiator, the internal mass flow, and the 
variation of the internal-flow characteristics with the addi- 
tion of heat. The effects of the addition of heat to the 
internal flow are also of importance because of the pos- 
sibility of obtaining appreciable thrust power fi'om the 
internal system. 

References 1 and 2, which are representative of the numer- 
ous theoretical analyses comprising most of the previous 
work in this field, were concerned mainly with evaluating 
the net force resulting from the flow through the engine- 
cooling system. In these earlier analyses it was usually not 
necessary to compute the internal-flow characteristics in 
great detail and it was possible to neglect, in part, the effects 
of compressibility (density changes). A more comprehen- 
sive general discussion of the cooling-system air flow was 
given in reference 3 but the details of the flow relations 
involved were not shown. 

In the present investigation a theoretical anlysis of the 
flow details at key stations in a typical system was made. 

of heating wei-e accounted 
Results of this analysis 
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Page 17^ coluim 1^ line 13 : The term at the beginning of the line 
should be • 



Page 3U, column 2^ line 31: The com-oression ratio should contain p 



P 

instead of p2; thus — = 1.226 • 
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The body used in this iiivcstiiralion was drsignod to con- 
form closely with the duct aiTaiigeinent assumed in the anal- 
ysis given in part I. Heat was added by a speciaUy con- 
structed radiator causing negligible blockage of the duct 
flow and insignificant friction loss. This manner of adding 
heat minimized the fire hazard m the wind tunnel and 
eliminated the combustion problem that woukl have existed 
if fuel were burned in the duct. The amount of licnt nddcd 
(160 kw max.) produced a thrust equal to about 40 jx^ccnl 
of the drag of the test body at M=0.75. This rate of heating, 
although lower than the rates obtainable by combustion of 
fuel, w^as large enough to permit accurate measurements of 
the thrust and other effects of heating to be obtained and 
thus to provide the basis for significant comparisons with tlic 
analytical investigation. 
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SYMBOLS 

area, square feet 

velocity of sound, feet })er second 
mass-flow coefficient ^ ^fy^ (^^") 

internal-drag coefficient (^^f^ 

internal-drag-coefficient incrcMuent (^^^^ 



drag coefficient of radiator ^^^^^^^(^^ ^ ^ 



radiator drag coefficient ( ^ 



at 



constant pressure 



(for air 0.24 



f(H4 



E. I 



specific heat 
Btu/lb/°F) 
!H^t drag due to internal flow, pounds 
internal-drag increment, pounds 
drag force due to friction within the radiator tubes 
pounds 

total drag force acting on the radiator, pounds 
duct diameter, inches 

maximum fuselage cross-section area, square 

(1.009 for model tested) 
acceleration of gravity (32.2 ft/sec^) 
heat added, Btu per second 
total pressure, pounds per square foot absolute 
total-pressure loss, pounds per square foot 
mechanical equivalent of heat (778 ft-lb/Btu) 
kinetic energy, foot-pounds 
over-all length of fuselage, inches 
distance from nose, inches 
Macli number (v/a) 
mass-flow rate, slugs per second 

pressure coeflficient ^ 

static pressure, pounds per square foot absolute 
static-pressure decrease, pounds per square* foot 
quantity of flow, cubic feet per second 

dynamic pressure, pounds per scjuarc foot (^^P^^^ 
gas constant, fcn^t per °F (for aii\ oo.:^) 
I'adius, inclu^s 



T 
T 
At 

Vo 

V 
X 

7 

er 

P 

Vd 



t(Mn])('rat urc, absohilc 
thrust, pounds 
temperature change, °F 
free-stream velocity, feet per second 
v(*lo('ity within duct, feet per second 
distance from leading edge of respective 
inches 

ratio of speeiflc hc^ats (for air, 1.40) 
heat-cycle efficiency (equation (31)) 
niechanical-effieiency factor (equation (33)) 
propulsive elliciency (€//e,u) ((npiation (32)) 
density, slugs per cubic foot 
diffuser efficiency 
compressibility factor 



sect ions 



14-— 4-^ + 
^ 4 ^40^1600 



Subscripts: 

0 to 5 stations in internal flow system shown in figure^ 1 
c test condition without heat 

d value at any station in duct 

/ friction component 

1 \o\v speed, incompressible-flow condition 

r condition across the radiator or heating device 

r2 within the radiator at the tube entrances 

rs within the i-adiator at the tube exits 

t condition W'ithin the radiator tubes 

A prime after a symbol indicates the condition without 
heat but with same mass flow as with heat. 

I. ANALYSIS OF INTERNAL FLOW SYSTEMS 

CALCULATION OF THE FLOW CHARACTERISTICS AT KEY STATIONS IN 
A TyPICAL INTERNAL FLOW SYSTEM 

In the design of efficient aircraft, the flow characteristics 
(static pressure, density, and velocity) must be computed 
at several key stations in each internal flow system. The 
objects of such calculations are to determine the state of the 
air entering the radiatoi's, the pressures available for cooling, 
tlu* flow changes across the heating device, the net drag or 
pr()pulsiv(^ force due to the internal flow, and the required 
sizes oi" inlet and outlet openings. It will be shown that 
under present-day operating conditions compressibility 
effects and the S(v*ondary efl'ects of heating cannot be n(^g- 
lected in the calculations. 

Methods will be presented for analyzing the internal flow- 
system that will account for the etfects of heating and 
compressibility and will yet be simple to use. Elinunation 
of much of the mathematical work will be effected through 
solution of the more cumbersome e(] nations in chart form. 
The uK^thods presented are generally applicable to any 
internal flow system, either with or without the addition of 
heat or mcu'lianical energy (internal blower), and for any 
rate of internal flow. 

One-dinuMisional flow and a uniform velocity distribution 
across the duct is assumed throughout the analysis. 

The theory of the process whereby the heat energy added 
to tlu^ internal flow^ is converted into mechanical energy is 
described under a subsecpient subheading of part I. 
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Initial conditions.— The initial quantities assumed to })e 
known b}^ the dosio-nor are as follows: 

(1) The rate at which heat is to be dissipated //, Blu 

per second 

(2) Tlie pressure drop lor the cold radiator {Aplq2)i 

obtained from low-speed tests 

(3) The mass-flow^ rate rec|uir(Ml by the system rn, slugs 

per second 

(4) The speed of flight Fo and tlie at niospluM-ic condi- 

tions />o, Po, and To 

(5) The duct areas available at all stations in the 

internal system 
Station 0. — In the method to be described, the internal 
i!()w is considered as a stream tube starting far ahead of the 
airplane where the velocity in the tube (relative to the air- 
plane) is Vq and the air conditions are those of the free stream 
(fig. 1). The area of the tube at this station is determined 
by tlu^ rc(juir(>(l rate of internal flow; that is, 



m 



(1) 



be 



The ar(*a ^lo is not(\'l to be a nn^asiire of the mass How and 
ars a simple relation to the flow coefficients, and ^ • 




Vo 


Vl ~ U„„r 


V2 


Po 


A 




Po 


A 




I 
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^0 


A, 





(a) 




£52Sn _ 




(b) P'or tubular radiiUor. 



(a) For internal-flow analysis. 

(c) For baffled air-cooled cylinder. 

Figure I— Station designation. 



If Q equals tlx^ Aolunie rate of flow at any station r/. 



The mass-flow coc^fficient C is 
C= 



m 

PqFV^ 



' F 



For isenlropic flow, the pressure, the (U^isity, and the 
temperature at any station in the duct for any rate of in- 
ternal flow can b(^ related to the known initial conditions 

through the area ratio A^^jAa and the Mach number Mo=-/; 

where a() = 49.0^To. Charts for obtaining these desired 
characteristics for known values of the area ratio and Mach 
number are given as figure 2 (pressure), figure 3 (density), 
and figure 4 (temperature). 

The mathematical l)asis for these curves was derivcul as 
follows: B}' B(Mnoulli's ('((nation for adiabatic compressible 
flow. 



1/2 
9 ' - 



7 Vo_v/ 



y Vd 



whe 



7—1 p(, 2 7—1 Pa 
From the requirement of flow continuity, 
PnAQVQ = pdAaVa 
Po^ 

Also, 

and for isentropic flow 



(2) 



Ptf/ V^tfy 



By substitntion of these relations into e([uation (2) in order 
to elinnnate the v(docity and density terms, the following 
rc^sult is obtaiiUMl: 



y— 

9 



(3) 



This equation is solved in figure 2 for pa/po in terms of 
Mo and Ao/Aa. The density and temperature ratios in the 
duct can be computed from the pressure-ratio values ob- 
tain(Ml from (equation (3). For isentropic flow^ 



Pd^ 

Po 



and 



7-1 
7 



The use of figures 3 and 4 permits the density and t(»mp(M"a- 
ture ratios to be obtained directly for the known vahu^s of 
Mq and Ao/Ad. 




Figure 2. — Variation of pressure ratio with area ratio for a given Macli number. DifTuser eiri(;iency= 1.0. 
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The use of the charts (figs. 2, 3, and 4) is not restricted to 
the Mach number A/o and the area ratio Ao/A,,y as indicated 
by the subscripts on the Hi>ur(»s. For examph*, if th(^ area 
ratio and Mach nunibcM- wcw taken as Ax/Af, and Mi, 
respectively, the charts would yield the ratios P]/pdf Pi.lpd, | 
and 7\/Tf,. Any other moi-e convenient area and corre- 
sponding Mach number might be similarly used in place of 
Aq and Mo. The charts may also be employed to determine 
the flow from a lai'ge area to a small area. Th(^ Mach num- 
ber used with the charts must always be \\\v Mach number 
corresponding to the smaller area. 

Station 1. — The entrance area Ai should be designed to 
allow the desired mass flow^ to enter the duct at an inlet- 
v(»locity ratio that will pcM-mit efTici(uit external flow as well 
as (»ffici(Mit intei'nal flow. lu the investigation reported in 
reCerence 4, velocity ratios of 0.3 to 0.6 were found to ful- 
fill these requirements for cii-cular openings of suitable 
sha()e located at the stagnation point of a body. Higher 
inlet-velocity ratios w^ill mako it diflicult to expand the in- 
ternal flow efficiently, and Iow(M' inl(4-v(4ocity ratios will 
generally make the external flow unsatisfactory for iiigh- 
speed airplaiK^s. If the inlet-velocity i\ is selected on the 
basis of tlH^s(* consid(M-at ions, th(» inh^t area is given by 



In most applications a sufficiently accurate value of ^li will 
be obtained when pi is taken as equal to po. 

After the entrance area and the ratio AMi have been 
determined, the ])ressure, the density, and the temperature 
at station 1 may be obtaiiUMl directly from the charts if the 
opcMiing is not located within a boundary layer. 

Station 2. — If the ducting losses between stations 1 and 2 
ai-e a jiegligible fraction of the free-stream dynamic pressure, 
the static pressure, the density, aiul the temperature at 
station 2 may be obtained directly from the charts (figs. 2, 
3, and 4) for the area ratio A0/A2. 

Under most conditions of flight, however, it will be neces- 
sary to account for tiu^ (lifrus(^r losses even in eflficient sys- 
tems. The characteristics of the diffuser are usually (h^fined 
bv a diffuscM- (^fliciency, 



or, pi'efei'ably. 




TJu' definitioji in terms of the static ])r(^ssures, as indicated, 
caji be used only in the reductioji of ijicomi)ressible-flow 
(lata because, at high s])ee(ls, the pressure lise is affected by 
factors other than the difl'user sha])e. The defijiition in 
terms of the total-pressure loss betwe(Mi stations 1 and 2, 
however, which yields the same value for at low specnls, 
remains essentially cojistant witli Mach luindxM- (s(h» refer- 
ence 5) and should therefore be used for all high-sjxuMl test 
data. If the velocity disti'ibut ion is i-easonably uniform, the 
total-pressure loss in the diffuser will e(|ual lli(^ decrease iji 
static pressure at statioji 2 below the ideal static j)r(»ssure for 
zero friction loss; that is, 

748(J24 — 17 2 



7 

Ap, = Mu qo) ( 1 - r),) = (7i [ 1 - (^4 ' ) ] ( 1 - ^./) 

The assumption of incompressible flow^ in computing the 
local dynamic-pressure ratio is justified in th(» usual case 
whei(^ the Mach number at the enti'ance is low and wluM'e 
A// 2 is small. 

The static-pressure ratio at station 2 is ther(dore 

V(i \Po/ ideal y>()L V^l,/J 

The ideal pressur(» ratio is obtained as Ixd'ore from figtu'c 2, 
and the diffuser efliciency 77,/ can be estimated from test data 
such as that given in references 5 and 0. Data on the diffuser 
effici(Micy for the model used in the present tests are given 
in reference* o. 

The tempcM'ature distribution across the duct at station 2 
is not uniform, the tempcM'at lu-e bcdng slightly higher near 
the walls where the friction (^fleets predominate*. A satis- 
facl oi-y value is given by 

To \Po/ ideal \P()/ ideal 

in whi(di the ideal pressure ratio of figure 2 is used. The 
corresponding value of the density is 

P2 = -||r = 0.000583 ^ 

The Mach nunilxM' of the flow in the duel at station 2 will 
b(* shown to be an impoi'tant i)ai'am(*ter and is giv(Mi by 

Flow changes within the radiator. -An exact deterniiiui- 
tion of the flow changes bcdween stations 2 and 3 must 
incliule calcidations of the flow details within the luxating 
d(»vice. The following analysis describes the flow through 
a tube or chaiuud of constant cross section in which lH\'it is 
added and in which fi'iction losses occur. 

If no appreciable (Miergy losses are assunu\l to occur at 
the tube (Mitrances, it is apparent that the flow changes 
between station 2 and station r, (at tlu^ tube* (uiti'ance) are 
det(M-niine(l by tlu* area ratio ArJAoand the "Mach number 
Mr,,, where Ar., is tlu* free ar(*a within the radiator. Figiu'cs 
2, 3, and 4 may therefore be used to obtain the desired con- 
ditions at station Vo by merely changing the nomenclature* 
on these charts. The charts require the use of tlu* 'Mach 
number A/,.,, which may Ix* comj)ut(Hl fi'om tlu* know 11 vahu* 
of Mo as follows: 

The vahu* of Mr., is (h^pcMident on p2/p,-^ and To/Tr.,, values 
which are to be obtained fioni the charts. A trial and error 
process must therefore be employed until corresponding 
values of P2/P,.,, and T.ITr., are obtained. Two or three 



J 
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trials usually sufRco. Tho pressure (li-o|) of the air cnteiinii; 
the tuhe is given hv 



(6) 



in whicli p2ipr2 is obtaiiu'd iliivctly i'roin figure^ 2, for the 
value of ArJA2 and Mr^ as just described. 

The pr(»ssure dro]) between stations r2 and (at the tube 
entrance and exit) is obtained by applying the momentum 
relation between these stations: 



where Df is the drag force due to friction witbin tlie radiator 
tubes and Ar^ = Ary The pressure (hop in terms of g^^ is, 
thus. 



^^2 



Qr^Ar 



(7) 



It is evident from equation (7) that the pressure drop in 
the tube is comprised of two components: First, a pressure 
drop associated with the frictional drag force and, second, a 
pressure drop due to the momentum increase in the tube. 
This second term, frequently neglected, becomes very im- 
poj'tant when the velocity in the tube is large. The drag 
term is presented in the form of a drag coefficient, denoted 
iiej-einafter as C^f, where 

The basic value of Cdj- correspontbng to the unlieated low- 
speed condition Cd^. can be determined fi-om pressure-(b'op 
measurements for this condition: 



This basic value of the drag coefficient depends on the 
geometry of the radiator tubes and on the Reynolds number. 
A consick^rable volume of pressure-drop data for radiator 
tubes is available and may be used to compute Cd^,. When 
heat is applied, the drag coefficient increases because of the 
increase in dynamic pressure along the tube. If this drag 
increase is assumed to be proportional to the mean dynamic 
prc^ssure in the tube. 



or 



?Kt)IK-S)] 



(8) 



In order to evaluate Cz>^ from the known l)asic value Cr,^^, 
the density ratia prjpr^ must be obtained. This relation is 
also required in evaluating ^Pt/qT2 (equation (7)). A method 
for computing prJPr^ exactly will now be described. 



The density ratio will be obtained by solving sinudtane- 
ously the momentum, energy, and continuity relations 
between the flows at stations /•2 and rz. Equation (7) shows 
the momentum relation. The energy equation, relating the 
total energy prr unit mass at the two stations, is 



- + 



7-1 



P:2,J_HJ!ii. 

PrJ m 2 



1 



(9) 



Since the areas at stations / j and /'3 are e(iual, 

Pr2^r2 = Pr/r., 

Solving the last three e(|uati()ns simultaneously for PtJpt^ 
gives 



CpgmT, 



(10) 



This relation has been evaluated for a wide range of values 
of heat input and radiator drag coefficients and the resiilts 
are given in figure 5 for use in design and performance work. 
Values of Prjpr^^ obtained from figure 5 with an estimated 
Cdj^ may be substituted in equation (8) to obtain a more 
exact value of Cd^- This process may then be repeated until 
the correct values of Prjpr^ and Cdj. are obtained. Two or 
three trials are usually sufficient. 

Several interesting conclusions may be drawn from a study 
of equation (10) or figure 5. When the Mach number Mr^ 
of the flow approaches zero, the value of pr.Jpr^ approaches 



1 + 



H 



c^gmTr 



or 1- 



the densitv ratio for heating at 



constant pressure. As Mr^ becomes greater, the pressure 
drop rapidly increases, with the result that the density ratio 
Prjpr.^ increases above the constant-pressure value 1 -) — - • 

This effect is very great for largc^ values of Mr,, such as are 
encountered with air-cooled engines. 

For each radiator drag coefficient a limiting value of M,.y 
is reached beyond which the solution of equation (10) be- 
comes imaginary. This limiting condition occurs when the 
flow velocity at the tube outlets reaches the speed of sound, 
that is, when Mr^=\.0. A curve connecting the limiting 
density ratios is shown in figure 5. When the drag coeffi- 
cients or the rate of heating is high, this critical condition 
is reached at low values of Mr^, owing to the large density 
decrease along the tube. No increase in mass flow can occur 
after this critical condition has been reached. It may be 
inferred from the general knowdedge of nozzles that shock 
waves will occur at the tube exits if the pressure drop is 
increased after sonic velocity has been reached. The energy 
corresponding to the increased pressure drop is dissipated in 
these shock waves and no increase in mass-flow rate occurs. 
The use of figure 5 as outlined will automatically indicate 
when this limiting condition is reached. 




2 



.3 



Mach numb en M^^ or 



(c) 



0.2. 



CpomiTr^ or 72) 
Figure 5— ContinucHl 
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With the density ratio and the drag coefficient evaluated 
by equations (8) and (10) as just described, the pressure 
drop between the entrance and exit stations of the tube is 
easily calculated from equation (7); that is, 



(11) 



Plow changes across the radiator. — The pressure change 
between stations 2 and 3 is the sum of the pressure drop due 
to acceleration of the air into the tubes plus the pressure 
drop within the tubes less the pressure recovery between the 
tube exits and station 3. The pressure-recovery term has 
yet to be evaluated. If station 3 is assumed to be far enough 
behind the radiator for the velocity distribution across the 
duct to be uniform, the momentum equation may be applied 
between stations 7-3 and 3 to give 



^.-^3=2,3(1-^;) 

or, as a satisfactory approximation, 

^^3-^3=-2^3(^^-l) (12) 

The over-all pressure drop across the radiator is the sum of 
equations (6), (11), and (12): 

Ap,=„-„.=„(l-^3)+,,.[c„,+<?5- 0 

The pressure drop is usually more conveniently expressed as a 
fraction of the dynamic pressure ^2, as follows: 



0.2 



(13) 



For the case of the tubular radiator, simple methods for 
evaluating ,2/2^2 and Pr2lp2 from figures 2 and 3, prjprz from 
figure 5, and Cof from equation (8) have been described. 
Equation (13) for the over-all pressure drop is exact except 
for the approximations made in simplifying the last term, 
which gives the pressure recovery between stations and 3. 
This term is a small fraction of the over-all pressure drop 
and cannot be computed precisely, in most cases, because 
th(i flow distribution rarely becomes uniform in the short 
duct lengths that are usually available for the pressure- 
recovery process. With air-cooled engines, for exam])le, the 
rear pressure station is always taken so close to the baffle 
exits that the pressure-recovery term is virtually zero and, 
hence, only the first two terms of equation (13) need be 
used. An example of the application of equation (13) to 
the case of a partly baffled air-cooled cylinder will be given 
in appendix A. 

Simplified method of computing radiator pressure drop. — 
The pressui-e chop across the heated radiator can be deter- 
mined by a less involved approximate method, accurate 
enough for many applications, which does not involve cal- 
culation of the details of the flow within the tubes or fins. 
The method consists of correcting the basic pressm-e di*op 



of the cold radiator to the design operating condition by a 
relation involving the density ratio P2/P3 across the radiator. 
The use of ps/ps instead of the density ratio within the tubes 
PtJptz eliminates the necessity of computing the tube flow 
and is justified by the fact that the two ratios are approx- 
imately equal. The same general procedure followed in 
deriving the exact equation for the pressure drop (equation 
(13)) will be employed in deriving the more simple approx- 
imate relation. If pa is assumed equal to p,2, the pressure 
drop between stations 2 and Vo is approximately 

Ap2-r2 = qr2-^2 = q2 [(^') " l] 
Within the tube (cf. equation (7) ), the pressure drop is 

or approximately 

The pressure rise between the tube exists and station 3 is 
(see equation (12) ) 

The sum of these three components is the desired ovei-all 
pressure drop ; that is, 

^--'".[i(-:;)]--(f:)X:-o 
+4(f:J-]-'.©(f:-0 <H, 

In order to use this reflation to compute the eft'ect of the 
density change aross the radiator, Ap;. must be evaluated 
from the known (or calculable) value of the over-all pressure 
drop of the cold radiator A/^r^. For the l^asic cold condition 
P2=P3, and equation (14) then becomes 

A,.„-Ay,.+4,[(^-J-l]-2„(^;-,) 
A,.=A,„-,.[(0'-,]+2,, (15) 



whence 



Thus Apf. can easily be computed from the over-all pressure 
drop for the cold radiator if the geometry of the radiator is 
known. The use of equations (14) and (15) to correct the 
pressure drop obtained for the cold condition (or for a test 
condition with heat) to the design operating condition has 
been found to yield results approaching in accuracy the 
results of equation (13) in cases where the Mach number of 
the flow in the tubes is not large. As previously discussed, 
the last term in both equations (14) and (15) should not be 
included in cases where station 3 is too close to the tube exits 
to permit the velocity distribution to become uniform. 
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In order to obtain the correct value of P2/P3 for use in 
equation (14), it is important that compressibihty effects, 
as well as heating effects, be allowed for. A simple method 
of obtaining P2/P3 from figure 5 will now be presented. For 
this method the drag coefficient of the radiator Cd, will be 
needed, and this value can be obtained from equation (14) as 
follows (cf. equation (11)): 



(7, 



^ Dr ^APr 



(16) 



where Apjq2 is given by equation(14). The flow in the constant- 
area duct between stations 2 and 3 is analogous to the 
flow in the tubes and exactly the same relations apply be- 
tw^een stations 2 and 3 as have been applied between stations 
Vo and r^. Equations (7), (9), (10), and (11) may thus be used 
to compute the density ratio P2/P3 and the pressure drop 
across the radiator by merely changing the notation to repre- 
sent the duct flow rather than the tube flow. Similarly, 
figure 5, which is a graphical representation of equation (10), 
may be used to obtain P2/P3 directly. The aforementioned 
changes in notation, applicable to equations (7) to (11) and 
to figure 5, are as follows: 

Replace Mr^ by M2 

Replace pjpr^ by P2/P3 

T) 1 friction r „ ^ radiator p p 

Replace ^ — by . ^ 01 Dy 0/>^ 



Figure 5 may be used with the revised notation to obtain 
a value for P2/P3. By use of a trial value of P2/P3, the drag 
coefficient Cd^ can be determined from equation (16). With 
this value of Co,, a new value of P2/P3 can be determined 
from figure 5. This process is repeated until corresponding 
values of Co, and P2/P3 are obtained. Equation (14) may be 
then evaluated for the radiator pressure drop. 

For design work, it is frequently necessary to compute the 
pressure-drop ratio Apjqo, This quantity can be written in 
the form 



q.0 g.2 \go/ 



which can be evaluated in terms of the flight Mach number 
and the design area ratio Aq/Ao in cases w^here the area 
ratio is small (say less than 0.2) and where the diffuser loss 
may be neglected. For these conditions, from equation (3) 
or figure 3, 

po=(^__j^ Y 

p, V1+0.20M0V 



\2.5 



Therefore, 



and 



2o P0VP2A/ V1+0.20M„7 VA/ 



APr 



22 



(17) 



Either the approximate value of Apr/qo given by equation 
(14) or the exact value from equation (13) may be used in 
equation (17). The quantities A0/A2 and Mq are design 
constants. 

It has been shown in figure 5 that the density ratio across 
a radiator or air-cooled engine becomes very large as the 
tube or duct Mach numbers increase, even with low rates 
of heating. Corresponding increases in friction and in the 
pressure drop due to momentum increase within the tubes 
have been indicated in equation (13) or (14). These adverse 
effects due to the density-ratio change should be emphasized 
because, under present-day operating conditions at high 
altitude, very large Mach numbers can be obtained within 
radiators and, in particular, within the baffles of air-cooled 
engines. Calculation of the required pressure drops will 
be seriously in error if no allowance is made for the eflects 
of density change. Equation (13) or (14) or equivalent 
relations must be used. High-speed test data are needed 
to indicate how closely the one-dimensional analysis of this 
report approximates the actual flow through an air-cooled 
engine. 

In order to illustrate the large pressure-drop increase 
that will be required as the operating altitude is increased, 
figure 6 has been prepared for a typical existing air-cooled 
engine. Figure 6(a) shows the variation of the pressure- 
drop ratio Apr/qo with altitude for the air-cooled engine. 



C>5 



1 I I r 

Sonic velocity 
within baffles 



Maximum available pressure 
drop without blower 




No allowance for density^ 
change across engine-.... 



ib) 



10 
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40 
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(a) Pressure-drop ratio in terms of djTiamic pressure ^PtIqo. 

(b) Pressure-drop ratio in terms of sea-level conditions. Afo=0.6 

Figure 6- Variation of pressure-drop ratio with altitude for a typical air-cooled engine. 
Engine conditions assumed are full-rated power; fuel-air ratio, 0.10; maximum head tem- 
perature, 450° F. 
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Two limiting conditions beyond which the desired cooHng 
conditions cannot l)e maintained are included. The first 
occurs when the maximum available pressure drop (assumed 
as l.Og'o with extended exit fhips) is reached. Flight at 
an altitude higher than the altitude defined by this limit 
would require the use of a blower. If a blower were used, 
the second limiting condition would occur at a still higher 
altitude, where sonic velocity would be attained within the 
baffles. A comparison is made in figure 6(b) of the increase 
in pressure drop with altitude as predicted by the methods 
of this report with the increase ])i'(»dicted by current methods 
in which no allowance is nuidc for the density change across 
the engine. Details of the calculation and the assumptions 
involved are given in appendix A. 

A study of the effects of increasing the fin area of the 
engine indicates that the cooling difricult ics sliowii on 
figure 6(a) could be delayed to altitudes higluM- than those 
indicated if the fin area could be increased. The same 
result could be obtained by the unsatisfactory exp(Mli(Mit of 
increasing either the fuel-air ratio or the head temperature. 
This fact undoubtedly explains the ability of some existing 
airplanes to operate at very high altitudes for limited 
periods of time. It is also evident that the cooling problem 
is alleviated by the reduction in power output that occurs 
if the critical altitude* of supcM'chai-ging is exceeded. 

Temperature change across radiator. — The temperature 
rise across the heating device is easily calculated after the 
pressure drop and density ratio po/pa have been obtained. 
The desired temperature-rise ratio MrlT2 is derived as 
follows: 

T^-T2 ___Mr__V,P2 ^ 
To To V2 P3 

wliere A2?r/g'2 is given i)y equation (13) or (14) and pa/ps is 
obtained from the chart (fig. 5) using the duct-flow notation. 

Station 3. — If either equation (13) or (14) is used to com- 
pute the ])ressure drop and equation (18) is used to compute 
the temperature change, the conditions at station 3 are 
gi ven by 

p., =0.000583 ^ 

Station 4.— The internal static pressure at station 4 
depends on the shape and location of the outlet. Efficient 



outlets of the type employed in the present tests are design (hI 
to make the streamlines of both the internal and external 
flows as nearly parallel as possible at the outlet; thus, the 
external outlet pressure becomes equal to the internal pres- 
sure. (See references 4 and 7.) The value of the outlet 
pressure for this type of outlet is found to be about the same 
as the pressure existing at that station on the body before 
the opening was added. Pressure-distribution data for 
faired bodies may therefore be used to estimate the outlet 
pressure for openings of the type described. If the walls of 
the outlet are converging rather than parallel, however, the 
outlet pressure will be considerably higher than the pressure 
on th(^ corresponding faired body. (See references 4 and 7.) 
With cusped outlets at the tail of the streamline body, as 
assumed for the present analysis, the outlet pressure is 
about 0.11(2(1 above stream static pressure, according to 
reference 5. If an estimate of the outlet pressure is made as 
suggested, the flow conditions at the outlet may be computed 
from 

P4 = 0.000583 

or 




The size of the opening required to permit passage of tlu^ 
design quantity of intei'iial flow can now ])e com])ut(Ml from 



where m is known, p4 has been givcMi previously, and V4 is 
given by 

=V(„^J-"::[-(;;:T1 

The assumption of isentropic flow between stations 3 and 4 
yields accurate results for unobstructed outlet passages. 
In ordinary cases the effects of friction are localized in a 
shallow boundary layer near the wall and do not affect the 
main outlet flow. In many practical arrangements, how- 
ever, the outlet passage is complicated by turns, tiers of 
guide vanes, and like arrangements, to which special con- 
sideration must be given. In such cases, it is suggested that 
the step-by-step process employed in determining the flow 
across the radiator be used. The drag coefficient of the 
obstruction can usually be estimated and figure 5 may be 
used to simplify the computations. 

A less exact but more simple method of accounting for the 
energy loss due to obstructions between stations 3 and 4 
assumes that the isentropic expansion takes place between 
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tho prossiin^ y>;^ — A/<3_.. and the outlet pressure, wIkmh^ 
Ahs-4 is the total-])ressur(^ loss oecurring- ])etween stations 3 
and 4; that is, as an approximation, 



P4 = P3 



A sample caU'idation usine this method was o-ivrn in 
j'(»ferenee 8. 

For usual installations, the maximum (hict velocity be- 
tween stations 3 and 4 will he attained at the outlet. It can 
he shown that the outlet Mach nmnher cannot exceed unity. 
This limitini^ condition cannot he reached ui ordinary cases, 
however, without the aid of an internal blower because the 
outlet Mach number will not exceed the flig-ht Mach mnubtM- 
exce])t in systems where the outlet is located in a resriou of 
high negative pressure. The addition of heat alon(^ will not 
appreciably change the outlet Mach iuiinb(M- bcM-ause, with 
the increase in outlet velocity due to the addition of heat, 
there is a corresponding increase in llie speed of sound. If 
th(^ critical vc^locity were reached in the outlet, the mass- 
llow rate could still be increased by an increase in outlet area. 

Station 5. — The velocity at station 5 is required for the 
])ui |)ose of computing the resultant force i\uv to the hiternal 
(low. As the static pressure is equal to y^i, by the general 
energy equation. 



-V(^,)"--"(';;)[-C:)'1 

The resultant force is obtaiiKMl by substitution of this value 
of into the monKMitum relation which, siiic(^ lh^ = po^ is 
simj^ly 



D = m\Vo-r5) 



A negative i-c^sult obviously 
forc(^ is a thrust. 



(li('at(>s that 



(21) 
resultant 



THE EFFECTS OF HEAT AND COMPRESSIBILITY ON THE TOTAL PRES- 
SURE AT THE OUTLET AND ON THE INTERNAL MASS FLOW 

Total pressure at the outlet. — In cases where the station- 
to-station nuahod is not emi)loyed, the total pressure at the 
outlet nuist be measured oi- calculatcMl in order to determine 
the net drag or thrust due to the internal flow. The station- 
to-station method of analysis just described, although essen- 
tial in calculating exactly the details of the internal flow, 
does not clearly indicate the over-all effects of heating and 
compressibihty on the outlet total pressure. A more direct 
method, which is useful in illustrating the effects in question 
and also in computing the mass-flow changes and the net 
thrust, will now be described. The total-pressure loss at the 
outlet is the sum of the loss across the radiator, the duct loss, 
and the pressure rise through the blower, if one is present; 
that is, 



\ (/o / oi er-all 



(22) 

\ ?0 / r \ 5o / duct \ (Zo / blower 



The total-pr(*ssui-e loss across the radiator is equal to the 
static-pressure droj) (equation (17)) less the dynamic-pressure 
rise (for the small vahies of A/2 usually encountered). 

\(loJr \qoJr qo \(Io/r (Jo\p3 / 

From (Kjuation (17), therefore, if no blower is used, 

+(-) 

V (Jo /duel 



(23) 



This (Hpiation is easily evaluated by use of figure 5 with th(» 
duct-flow notation to obtain the necessary values of po/pg. 
The duct loss should include an estimate of the skin-friction 
losses plus the difTus(*r loss, which, from equation (.">), is 
given by 

\qo/di/ruser goL \^d/ J 

For a constant mass-flow ratio A0/A2, the radiator drag 
increases as a result of the incrc^ase in density ratio P2/P3 due 
to heating. As shown by equation (23), the total-pressure 
loss will correspondingly increase. Aside from this well- 
know^n effect, howt^vcM', the addition of heat will cause httle 
change in tlu^ outlet total pressure. From the development 
leading to equation (23), it is seen that, for a given mass 
flow and i-adiator drag, the outlet total pressure will de- 
crease as heat is added by the amount q^—qo. Ow'mg to the 
low duct velocities, the quantity {q^ — q2)lqo will be very small 
in usual cases. This result is of considerable interest in that 
measurement of total pressure in the wake will give no 
direct indication of the propulsive effect of heating. 

When the duct velocities are large, however, and the 
corresponding value of 53—^2 is appreciable, the foregoing 
d(»v(^lopm(Mit indicates that the eft'ect of heat will cause an 
appreciable loss in outlet total pressure over and above the 
increascnl radiator drag loss. 

For a duct system with moderate intt^nal losses, no ap- 
preciable drag change of surfaces exposed to the outlet flow 
is to be expected because of the negligibly small chang(» in 
outlet dynamic pressure with heating. The Reynolds num- 
ber of the outlet flow^ decreases, however, and some scale 
effect might therefore be expected. No change in the out- 
let Mach number occurs because the velocity increase is 
offset by an increase of the speed of sound in the heated 
outlet flow. 

Variation of mass flow. -The internal nniss-flow coeflicient 



m 



-pPTr- or the (low ratio is a constant for a given 

PqFVq Po v 0 

system, except for the effects of compressibility and heating. 
(See reference 5.) In order to show these effects, the ratio 



will be written in terms of Mo and 

Povo 



H 



From th( 
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general ciu'i-gy equation, 



PiV4 

9 



whence, 



where 



and 



PoVo Y 7— 1 Po qo [KpJ J 



(24) 



The density ratio p^/po in equation (24) can be written in 
terms of the lieat added and A/o. 



P4 _ P5 P4 _ n (Pj\ '^ ^^(1+ 

Po Po Po T-o \po/ T^x 2 / 

If the internal losses are assumed to be moderate and the 
pressure at the radiator to approach stagnation pressure, it 
can be shown from the energy relations that 7 5 is related to 
the heat added as follows: 



To 



1 + 



H 



wlu'nco 



P4_ (I+O.7P4M0')''"* 
Po , , H 



1 



(25) 



7.74TOro(l+0.20iW) 



When the indicated substitutions for total pressure are made, equation (24) becomes 

/ 



PoVo \ po-v/ - 



7 (2±yRMl 
-iV 7Mo^ 



) 



(2 + 7P4MoO( 



1- 



1 



J 1 2+7P4M0 ' 



0^ 



(26) 



where p4/po is given by equation (25). Tiis])ection of this 
result shows that, for fixed values of il/o i^i^d Ah/qo, the mass 
flow will decrease as heat is applied in proportion to p4^^~. 
The outlet area, therefore, must be increased as heat is addend 
to the flow if a fixed mass-flow rate is to be maintained. 
Equation (26) has been solved for a range of values of 
Ah/Qo and Mq in order to show the effect of Mach number and 
to facilitate the use of this relation in design work. The 
solution is shown on figure 7 for an outlet pressiu-e equal to 
free-stream static pressure ^4=0. The Mach number effect 
on the flow ratio is seen to be appreciable. Equation (26) 
or figure 7 will yield values of p^r^ as accurate as the results 
of the step-by-step process for most applications, provided 
that the internal losses are not excessive and that the heat 
is added at a duct pressure approaching stagnation pressure. 
The values of Ah/qo aiid p^/po required for solution of 
equation (26) or for use with figure 7 are obtained from 
equations (23) and (25), respectively. 

If the flow is considei-ed to be incompressible except for 
the effect of heating, equation (26) can be reduced to 



PoVo Vpo/ V go / 



(27) 



as Mq approaches zero. This relation was given in reference 7. 
The curve on figure 7 for Mo=0 is a plot of equation (27) 
for the outlet pressure coefficient Pi=0, Because of its 
simplicity and also because equation (26) tends to become 
inaccurate at very low Mach numbers, equation (27) should 
be used for Mo<0.2 in preference to equation (26) in cases 
where figure 7 is not used. Figure 7 applies for only one 
value of the outlet pressure coefficient Pi = 0, It may be 
used for values of F4 other than 0 by correcting the result 
obtained from the figure to the required value of P4 according 
to equation (27). 

Direct calculation of the outlet area. — Calculation of the 
size of outlet opening required to pass the design mass flow 
is one of the most frequently perform.ed operations in the 



design of an internal flow system. If the step-by-step 
process has been carried out, the outlet area may be calcu- 
lated from the values of p4 and V4 as already shown (equation 
(19)). In general, however, the computation of A4 without 
the detailed intermediate steps is desirable. Equation (25), 
with equation (26) or (27), or figure 7 can be used to calculate 

— y which peril 
Po V 0 

. m m. 



the outlet area to be oblaiiied directly: 



P4^'4 



or 



\Po V 0/ 



Pol 0 



/ PaVa \ 



(28) 



Effect of an operating propeller. -If the inlet opening is 
located in the slipstream, I lie velocity, the pressure, and the 
density in the slipstream should be taken as the initial con- 
ditions and the same procedure followed as for the condition 
without propeller in arriving at the conditions hi the duct. 
In cases where the inlet is immediately behind the propeller, 
it will be found convenient to take the flow characteristics 
in the inlet opening as the initial conditions rather than the 
conditions ahead of the propeller. The area-ratio charts 
(figs. 2, 3, and 4) may be used after substitution of and 
Ml for Aq and Mq. 

THE CONVERSION OF HEAT ENERGY TO THRUST POWER 

The effects of the addition of heat in increasing the 
radiator drag and in reducing the internal mass flow have 
been discussed. Both of these eftects cause changes in the 
drag due to the internal losses, but neither effect results 
from the recovery of mechanical energy from the heat added. 
In the study of the propulsive eflect due to heating, there- 
fore, it is desirable to eliminate these unrelated drag changes. 
This elimination can be accomplished by assuming, for 
comparison with the heated condition, a hypothetical con- 
dition without heat in which both the radiator drag and the 
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Total-pressure /oss, Ah/q^ 

Figure 7.— Cui'ves for determining the outlet velocity rati( 
(Solution of equation (26) for P4=0.) 



.7 



mass-flow rate have the same values as in the heated condi- 
tion. This assumption not only simplifies the analysis of 
the recovery process but also gi'eatly facilitates calculation 
of the net propulsive force, as will presently be shown. In 
the following analysis this assumed unheated condition will 
])e indicated by the use of primes with the symbols. 

The net drag due to the flow tlu-ough an internal system is 
equal to the decrease in the momentum of the internal flow, 
provided that the momentum is measm-ed at stations where 
the static pressures are equal to free-stream pressure. 

From equation (21), 

If the system produces a net thrust, this thrust is given by 
th(» inci-(»ase in momentum of the internal flow: 

T=-D=m{v--V,) 

I'he (h'ag decrease (or thrust increase) resulting from the 
rcH'overy of mechanical energy from the addition of heat is 



Ar= (Thrust with heat) — (Thrust without heat but 
with same mass flow and same radiator drag as 
with heat) 



or 



^T=m{v^-v^') 



(29) 



The increase in kinetic energy of the internal flow due to the 
recovery process is thus 

A(K E^^'^iv.'-vn (30) 

If the efficiency of the process is defined as 

K. E. recovered 
^''^ Total heat added 

Then, from ec^uation (30), 



2JH 



(31) 



The over-all efficiency of conversion of heat cjiergy uito 
thrust power is the ratio 



Useful thrust pow er _ATVo_ m{vs—V5)Vo 



Heat add 



JH 



JH 



(32) 
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This oqiiatioji can bo writtoji in torms of tin* heat-cycle 
efficiency (equation (31)) as follows: 



(33) 



The (luantity e^f represents the efficiency witli wliicli tlie 
kinetic energy recovered from, the heat input is converted 
into thrust power. This efficiency is analogous to the ideal 
efficiency of a propeller and, in the same way, accounts for 
the energy left in the wake when the wake velocity is different 
from Vq. It will be found convenient in the succeeding 
analysis to express in terms of the internal-drag 
coefficients with and without heat. From equations (21) and 
(33), expressing the (hag and mass flow in coefficient form, 



1 



1 (Cr>~\~Cn) 



(34) 



Efficiency of the heat-recovery cycle. — If the duct velocity 
at the radiator section is very small in com.parison with the 
flight velocity, the static pressure at which heat is added 
will approach stagnation pressure and the pressure change 
across the radiator will approach zero. Under these condi- 
tions the heat cycle consists of an adiabatic compression 
from the stream static pressure to stagnation pressure, an 
addition of heat at constant pressure {p2^Pr^Ps)) and an 
adiabatic (expansion back to 2>o. The continual admission of 
air at pressure 2^0 ^ind density po is equivalent to cooling the 
outlet flow back to po at constant pressure po and this step 
completes the cycle. This process corresponds to an ideal 
heat engine cycle, the efficiency of which is, from elementary 
therm.odynamic relations, 

1—7 

_ Kinetic energy recovered _ f P2\ ^ /or\ 
Total heat achh^l ~ ^ ~ W ^ 

The assumptions on which this formula is based are justified 
in the case of the internal flow system only when the duct 
velocities are very small by comparison with the flight 
speed and when the friction losses in the radiator or combus- 
tion chamber are negligible. Operating conditions exist, 
however, in which relatively high duct velocities occur, and 
an expression for the efficiency applicable to such cases will 
now be derived. The heat-cycle efficiency will be shown 
to depend not only on the compression ratio but also on the 
Mach number of the duct flow and on th(^ rate of heatinfr. 



From the energy equation (cf. equation (9)), 

2 m ' 
By equation (31), therefore, 



-Jgc,{T,-T-J) 



(36) 



The values of T5 and T^^ will depend on the outlet pressure 
and on the conditions at station 3, which in the general case 
will depend on M2 and on the radiator or combustion-chamber 



losses. The previous analysis has indicated that when a 
radiator of conventional design is present it is impossible^ 
to obtain high Mach numbers in the duct at station 2 Ix*- 
cause the flow will choke within the radiator at tlu^ tube 
exits for relatively low values of Mo. Since it is d(»sii'(nl in 
the present analysis to show the effects of relativc^ly large 
increases in Mo on the (efficiency of the recovery cycle, it 
will be assunKvl that no radiator obst/-ucts the duct. The 
addition of heat will be assumed to take place by combustion 
of fuel in the duct between stations 2 and 3 by means of a 
burner having negligible friction losses. This ideal arrange- 
ment is ])eing approached in the ram-jet types of propulsion 
system. The effects of the increase in mass flow due to 
addition of fuel are of secondary importance and will be 
neglected. From these assiunptions T^^ can be related to 
To as follows: 



The pressure 7)3 is less than po by the value of the monuMitum 
increase due to the addition of heat; that is, inasmuch as 

^^2 — ^3 



PZ=P2 — P2V2- 



P3 
Pi 



= 1-7^2^ 



(2-) 
(:-') 



Also, fj om the equation of state, 



P2PZ 



P3 P2 



and, substituting for p^lpz^ 

T,= Tj'\\-y\l 
P3 L 



Now, 



=^'i{ff[-M%-')Y 

Tlu' licat-cyde efficiency is, therefore, from equation (36) 



\cpgmTj 



(1-0 



-1 (37) 



This equation for en is plotted in figure 8 as a function of 
compression ratio for various values of the heat input factor 
and M2. The term pa/ps in equation (37) is a function only 
of these variables and was computed from equation (10) for 
(7^ =0, using notation appropriate to the duct rather than 
the radiator tube. 



Flight Moch number-, Flicjht Mach numbGr, 

0 l.O 1.2 U 1.6 1.7 O 1.0 1.2 1.4 1.6 1.7 




Figure 8.— Variation of hoat-cycle einciency with pressure ratio and flight Mach number for various internal Mach numbers at station 2. C£>^=0. 
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Figure 8.— Concluded 
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The Mach nuinber at the entrance to the combustion 
cliamber, M2, is shown in figure 8 to have a large effect on the 
heat-cycle efficiency. When M2 is very small, approaching 
zero in value, the heat-cycle efficiency approaches the value 
given by equation (35), regardless of the heating rate. As 
M2 increases, however, the value of decreases because of 
the pressure drop due to acceleration of the heated air in the 
combustion chamber. For a given value of il/2, increasing 
the rate of heating causes decreases in €^ for the same reason. 
The limiting or ''choking" values of M2 are also shown in 
figure 8. If no heat were added, the highest possible value of 
M2 would be 1.0. As heat is added, the choking value of 
M2 decreases because sonic velocity is reached at the end of 
the combustion chamber {M^=l.O) for values of M2 less 
than 1.0. At the high rates of heating (for example, at 

c gmT ^^'^^ which corresponds to a temperature rise of the 

order of 2500° F in the combustion chamber), the choking 
value of il/2 becomes very small, only 0.196 for the example 
cited. This efl'ect is extremely important because it dictates 
the minimum possible area expansion ^2Mo for a given 
heating rate. 

The flight Mach numbers required to produce the com- 
pression ratios of figure 8 (assuming ideal stagnation pres- 
sure is obtained at station 2) are noted on the figure. At 
the higher supersonic speeds, the ideal thermal elhciencies 
for low values of M2 approach those obtainable in internal- 
com])ustion engines. This result would be expected in view 
of the fact that the compression ratios are of comparable 
magnitude. 

Because en increases with compression ratio, the inlet 
ducting and diffuser losses should be as low as possible. If an 
internal blower is not used, high pressures in the heating 
section of the duct can be obtained only when the area of this 
section is greater than the initial area Aq of the internal-flow 
tube. (See fig. 2.) This location of the heating device in a 
high-pressure, low-velocity section of the duct is also ideal 
from the standpoint of minimizing the losses in the lunitrng 
device. 

The friction loss within the heating device, A/i„ has the 
undesirable efl'ect of reducing the compression ratio. For 
reasons of simplification this loss was neglected in the evalua- 
tion of e//, equation (37), figure 8. If the principal part of 
this loss occurs upstream from the region in which heat is 
added, it may be accounted for by using ^2— A/ir instead of ^^2 
as the effective pressure at the entrance to the heating region. 
In the case in which the friction loss is intinuitely involved 
in the heating process as in a tubular radiator, the actual 
pressure drop during heating will be greater than assumed 
in equation (37). As an approximation in this case, it may 

be assumed that the effective pressure is ^2 ^ at the 

entrance to the heating chamber, and this value may then 
be used in equation (37) or figure 8 to obtain e/^. In l)oth 
of these examples, the value of A/i, should include estimates 
of the heating and compressibility effects on the friction loss 



if low-speed data for unheated devices are used in the 
calculation. 

Calculation of the net internal drag (or thrust). — In per- 
formance studies, it is desirable to compute the resultant 
drag, or thrust, due to the heated internal flow without going 
through the step-by-step duct calculations previously dis- 
cussed. The method to be described is exact within the 
limits of the assumption that the outlet flow is uniform. 
Very few steps and only simple slide-rule operations are 
involved. 

The net force due to the internal flow (equation f21)) is 
Z}=m(Fo-'i^5)=m(Fo-^50-m(t;5-?^50 (38) 

Equation (29) shows the significance of the two terms of 
equation (38) to be: 

(Net force due to internal flow)=(Internal drag 
without heat) — (Drag reduction due to heating) 

Each of these components can be computed without dlfii- 
culty, and the separation of the two is the basis for tlu* 
present method. 

If the flow were incompressible, the (h-ag corresponding 
to the unheated condition could be very easily computed 
from 

Z>/ = m(F„-0=m[r„-(TV-^^|)"'] 

or 

where Milq^ is obtained from equation (23) for values of 
mass-flow coefficient and radiator drag corr(^sponding to the 
heated-flow condition. For compressible flow, unfortun- 
ately, the drag calculation is laborious. A relation exists, 
however, between the ti'ue value of Cd for compressible flow 
and the value computed from the simple incompressible-flow 
formula, equation (39). From equations (21) and (39), 




This equation involves only tiie total-pressure loss A/i/go 
and the velocity ratio z^s/T^o. which can be defined in terms of 
A/i/(/o and from equations presented in appendix B. A 
plot of this relation is given in figure 9 and may be used to 
obtain the true value of the internal drag for the unheated 
condition using the value of the low-speed drag coefficient 
computed from equation (39); that is, 

the quantity being taken from figure 9. 
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Figure 9.— Ratio of tlio interiial-draj: cooniciont for compressible flow to the value computed for incompressible flow as a function of Mach number. 



Tlie drag for the iinliGated condition having been already 
obtained, only the drag reduction due to heating (second term 
of equation (38)) remains to be computed. From equa- 
tions (32) and (33) the drag reduction, or thrust, due to 
heating is 



or 



iloFVo 



(40) 



In summary, the net or resultant force (taken as positive 
in the drag direction) is given by the following formula: 

0. = C,,--AC.-2C[1 -(l -0'](g0- ^..„ (41) 

in which the parameters involved are computed as follows: 
Ah/qo equation (23) 
{CdICd.Y from figure 9 

en from figure 8 or equation (35) or (37) 

€m from equations (34) and (41) by trial 

The separation of the propulsive effect of heating from 
the internal drag due to frictional effects is useful in several 
other connections. In the correction of wind-tunnel (hita 
obtained without heat, for example, the right-hand term of 
equation (41) may be used to evaluate the (h-ag lechiction 
due to heating. It is obvious that model t(*sts should be 
made with values of Ah/qo aiid C (•()ri('sj)()iHrnig to the 
heated condition. A further application of the separation 
principle will be made in part II, where the thrust force due 



to the added heat is separated from the net force due to 
the internal flow. The results analyzed on this basis are 
significant as applied to any internal system. 

II. HIGH-SPEED WIND-TUNNEL TESTS OF A MODEL 
RAM-JET PROPULSION SYSTEM 

APPARATUS AND METHODS 

The Langley 8-foot high-speed wind tunnel, in which 
these tests were carried out, is a closed-throat, circular- 
section, single-return tunnel with airspeed continuously 
controllable for specMls from about 75 to 600 miles prv hour. 
The turbulence of the air stream is unusually low but is 
somewhat higher than in free air. 

The mo(l(4 tested consisted of a 13.G-inch-diameter body 
of revolution housing a 160-kilowatt heating device and 
mounted on a 9-percent-thick, 20-inch-chord, low-drag 
airfoil that spanned the jet (fig. 10). The air inlet was at 
the nose of the body with an expanding duct section back to 
the radiator. The air outlet was at the tail. The general 
model layout is shown in figure 11 with the ordinates listed 
in tables I and II. 

Th(* external shape of th(^ wing-fuselage combination was 
designed to have a critical speed in excess of the maximum 
test speed. In order to meet this criterion, inlet and outlet 
openings designed according to reference 4 were utihzed, 
and a O-percent-thick 20-inch-chord wing of airfoil section 
NACA 66-009 was required. The inlet- and outlet-opening 
shapes of reference 4 had the further advantage of permitting 
wide variations in the rate of internal flow and in the inlet- 
and outlet-velocity ratios without causing appreciable 
(>xternal-drag changes. 
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Figure 10.— Model assr;mbly. 
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Airfo// oncl/nate5 NACA66-009 
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.209 


If .00 


.f^\ 


.50 


.282 


12.00 


.642 


1.00 


.389 


13.00 


.762 


1.50 


.469 


14.00 


.69/ 


2.00 


.535 


15.00 


.568 


3.00 


.636 


16.00 


.444 


4.00 


.7/6 


17. OO 


.319 


5.00 


.762 


18.00 


.185 


6.00 


.632 


19.00 


.074 


1.00 


.366 


20. OO 


0 


Leading- edge radium =0.1 12 



/3.5 -\i5.o- 



Fuse I- J je ordina^e5 given T 
in tables 1 and E 
All dimen5ioii5 in inched 




13.5 in. Station 



/3.65 —J 



Part/al -annular outlet , tail E 



- 20.0 - 



Thermocou/Dle rake 
Pressure rake 
Portial-annular outlet, tail E 
Streamline tail 

^ Tail rake 




Fkjurk 11.— Model details, 



TABLE I.— STREAMLIXK-1U)DV ORDINATKS IN INCHES 



Streamline nose 


Center section 


Streamline tail 


X 


R 


X 


R 






0 


0 


0 


6. 78 




0 


6. 49 


.73 


1.29 


1.5 


6. 80 




3.0 


6.25 


1.46 


1.95 


4.5 


6.80 




6.0 


5. 92 


2. 92 


2. 93 


10.5 


6.80 




9.0 


5. 52 


4. 37 


3.64 


13.5 


6.80 




12.0 


5.06 


5.83 


4.21 


15.5 


6. 78 




15.0 


4. 52 


8. 75 


5.03 


18.5 ' 


6.71 




18.0 


3. 94 


11.66 


5.64 


21.5 


6. 58 




21.0 


3. 35 


14. 58 


6. 10 


23.0 


6. 49 




24.0 


2. 75 


17.49 


6. 44 








27.0 


2. 14 


20.41 


6. 67 








30. 0 


1.53 


23. 00 


6. 78 








33. 0 


.92 












37.5 


0 



Three ciisped-type openings at the tail (figs. 10 and 11) 
were employed to vary the internal mass-flow coefficient 
C from about 0.07 for the largest size (tail A) to about 0.025 
for the smallest size (tail C). Tail D was similar in size to 
the intermediate cusped outlet (tail B) but had straight, 
slightly converging walls. Tail E was a partial-annular 
outlet located at the 56-inch station on the body and d(^- 
signed according to reference 7. Tails B, D, and E had 
approximately the same flow coefficient, namely, 0.052. 



TABLE II.— DUCTED-HODV ORDIXATES IN INCHES 



Diffuser 



0 

. 10 
.25 
.50 
.75 
1.00 
1.75 
7. 75 
14. 75 
16.75 
18. 75 
20. 75 
22. 75 
25.00 



2. 34 
2.21 
2. 21 
2. 21 
2. 22 
2.23 

2. 28 
2.69 
3. 17 

3. 34 
3.59 
3. 98 
4.53 
5.59 



Nose radius: 
0.10 



Ducted body 
nose ordinates 



0 

.025 
.05 
. 10 
.25 
.50 
.75 
1.00 
1.75 
2. 75 
3.75 
4.75 
6. 75 
8.75 
10. 75 
12. 75 
14. 75 
16. 75 

18. 75 

19. 50 



2. 34 
2.44 
2. 50 
2.58 

2. 75 
2. 98 
3. 16 

3. 33 
3. 73 
4.15 
4.49 
4.79 
5.29 
5.71 
6. 04 
6. 30 
6. 51 
6. 66 
6. 76 
6. 78 



Tail 



0 

3.00 
6.00 
9.00 
12.00 
15.00 
18.00 
21.00 
24.00 
27.00 
28. 75 



Outlet 
diam. 



Outlet 
area 



6. 49 
6.25 
5. 92 
5. 47 
4.89 
4.24 
3.57 
2. 95 
2.42 
2.07 
1.985 



11.07 



6. 49 
6.25 
5.91 
5.42 
4.80 
4.09 
3. 37 
2. 70 
2. 17 
1.82 
1.705 



6. 49 
6. 25 
5.88 
5. 37 
4.71 
3. 96 
3. 18 
2.46 
1.87 
1.48 
1.365 



!.68 



The internal system was designed to house the radiator 
and survey equipment with a minimum of duct losses. A 
more detailed description of the model and design consider- 
ations is given in reference 5. 



HKA'I' AND COMPHKSSlHII.rrV EFFECTS IN INTEHNAL FLOW SYSTEMS 



29 




I'K.L lu: 12. -Uadiat.or details aiul inslallalion. 



RADIATOR 



Tlu* hoat input to the intornal air How was provided by a 
l()0-kilowatt, 220-volt, throo-phase electric resistance radiator 
designed for a maximum operating temperature of 1500° F. 
The radiator consisted of six elements of V/i inch X 0.020 
Nichrome V furnace strip (radiating surface, 32 sq ft) woven 
on a framework of reinforced asbestos millboard supports 
(figs. 12 and 13). Radiator temperatuiv was measured by 
four chromel-alumel tluM-mocoupl(^s welded to the elements. 

Power input to the radiator was controlled by eight trans- 
former voltage taps spaced to provide approximately uni- 
form power incremcMits. Total power input was measured 
b}^ a standard three-phas(^ wattmeter. 



SURVEY EQUIPMENT 



Pressure surveys. — Pressures ahead of and behind the 
radiator wei-e uK^asuied, respectively, by an 8-tube and an 
11-tube static- and total-pressur(» rake spanning the duct. 
(See figs. 11 and 12 for installation details.) All pressure 
tubes were connc^'ted to a standard manometer and the 
pressure readings were photographically recorded. A 52- 
tube removable survey rake of static- and total-pressure 
tubes was installed at station 4 for measurement of the 
internal and external pressures at the tail outlet. The i-ake 
tubes were supported by two 2}^inch-chord blades (NACA 
section 16-010) set at right angles to each other and extend- 
ing beyond the outlet walls (fig. 14). The tail rake was sup- 
ported by a l}^-inch outside-diameter steel tube housing the 
pressure leads and extending through the core of the radiator 
to the wing duct (fig. 11). 
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^^<:f'sfo/- supports--. 

JaU-rahe tube-. 




I''i(;r !{E 13. Kadi'itor clciiiciii . 




FiGUKE 14.— Tail-rake installation. 




Figure lo.— Thormocouplo rakr. 



Temperature measurements. — All temperatures were meas- 
ured l)v ('alii)rat(Ml thormocouph^s counoctc^d to a sensitive 
potentiometer. Air tenij)ei'ature ahead of lh(* radiator was 
measured by two shielded thermocouples located in the front 
th(»rmocoupIe rake. (See fi^:. 15 for rake details.) 

Shieldin;:^ of the junction prevented false air-temperature 
readings due to radiation effects. Air-temperature rise across 
th(^ radiator was measured dircH'tly l)y two thermocouple 
rakes containing six junctions in separate com])artnients 
(fig. 15). The tliermocouples in the fiont rake were tlie 
cold junctions and the corivsponding thernio(OU])l(»s in tlu^ 
rejir rake were the hot junctions. 

Temporary apparatus. — During preliminary calibration 
runs, additioiml apparatus was installed in the model to sup- 
plem(Mit the fixed e(juij)ment. A 14-tube total-pressure rak(» 
spanning the duct was mounted behind the radiator to aid 
in determining tlu^ duct velocity profile. From these data 
th(» energy addcul to the air stream was determined b}' an 
integration of the mass flow and the air-temperature rise. 

The heat lost through the walls of the tail was computcul 
from measurements of tlu^ wall temperatures and local skin- 
friction eoeflicients. A standard 5-tube static- and total- 
])r(»ssur(^ ^ 'mouse" was used to survey the boinidaiy layer of 
tail B in conjunction with local wall-temperature measure- 
UKMits ])y th(M-ni()('oii])l(^s. 

TESTS 

Drag tests with each outlet \v(M'e carried through a Mach 
numb(4- I'ange of {o 0.7.") i ;i ppiox. 200 to ")()() m])h in the 
I 8-ft higli-s])eed tunnel). tail rake (fig. 14) was removed 
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I. 
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□ 

X 
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A 
A 
A 
A 
B 
B 
B 
C 
C 
D 



Heat input, H/A^ , Blu/sq ft/sec 
A-r =-6St Sq ft 



9/ 
105 

142 
/66 
65 
104- 
128 
46 
96 
85 

B, baff/e 103 
plate 



r 



4t 



/.2 



Kir.rHE U). Varialion of hcat-cycic rllicicncy with pressure r 



/.3 

ratio , p^/po 

lit io Tor various heat inputs. 



Data from pressures at outlet. 



/.5 



(Im-inir \ lorcr t(»sts. Internal pressures at stations 2 and 3 
were recorded simultaneously with the force data, as was 
the power input to the ladialoi'. The range of power input 
was from 0 to HiO kilowatts. 

Tests of the model with the tail outlets A, B, C, and D 
were made with the tail rake ffig. 14) in place and were similar 
in range to the force tests. The object of these runs was to 
furnish data for calculations of the Ikn-U encM'gy recovery 
independent of forc(» mc^asuriMuents and also to evaluate 
external-drag changers under various heating conditions. A 
special run was made with tail B with an int(M'nal resistance 
plate that ])roduce(l a total-pressure drop of 0.27go. 



Supplementaiy I'uns to j)ermit an estimation of th(» heat 
lost througii the walls included measurements of wall 
temi)eratui'e and boundary-layer velocity profiles. 

All tests w(M'e made at zero angle* of attack. 

RESULTS 

The methods (employed to compute the lu^at-cycle 
(^dieiency and tlu* duct-flow chai'acteristics from tin* test 
data are described in aj)j)(Mi(li.\ B. 

Figure 16 shows the \ ariation of heat-cycle efliciency with 
compression ratio as determined from the wake-survey data. 
The experinuMital results are shown for all the tests and there- 
fore cover a wide I'ange of mass-flow rates, heat input, Mach 
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llr 



lO- 



Heut input, H/Ar , Bfujsq ff/i 

Ar~ =68lsq ft 




1.2 1.3 
Pnessune ratio jp^ 

FiGTRE 17.— Varial ion of heal -cycle clliciciicy with i)r('ssure ratio for various heat inputs. 



1.5 



Data from force nicasurfmonts. 



iiuin])or, and diictiiuj: losses. The lioat-cych* ofricioncv iiidi- 
caU'd by tlio forc(»-t(»sl data is sliovvii on fio-urc* 17 for the 
same ran^^e of varinhlcs as presented for the pressure 
measureineuts. 

Comparisons of tlu^ (low eharaetiM'istics measurtnl at key 
stations in the duet with the ehai'aeteristics computed by the 
methods of part 1 are shown in figui'e 18 Table 111 shows a 
comparison of the test values of total-])ressure loss with the 
corresponding values computed from th(^ theory. Table IV 
shows a comparison of experimental and theoretical mass- 
flow ratios. Figure 19 com]3ar(»s the (^xpiM'inuMital variation 
of mass-flow coefficient with outlet density ratio with the 
variations computed from the sim])lifi(Hl theory (ecpiation 
(27)) presented in part I. 

DISCUSSION 

Eflficiency of heat-recovery process. — The maximum Mach 
numbers of the duct flow at the j'adiator section did not 
exceed 0.10. K(^ference to figure 8 or to equation 0^7) will 



show that, for this ]ow value of .\/2, the heat-cycle efficiency 
depends only on the compression ratio yilv^'^ ^^^^^ is, 



6//=l- 



TABLK 111.— COMPARISON OF EXPERIMENTAL AND 
THEORETICAL VALUES OF TOTAL-PRESSURE LOSS 

[Tail B, baffle plate] 



M 


Total-pressure loss, M/go 


//= 


=0 


//«90 Btu/sec 


Equation 
(23) 


Exix^rimon- 
tal 


Equation 
(23) 


Experimen- 
tal 


0. 263 
.347 


0. 299 
.299 


0.301 
.298 










. 403 


.298 


.296 


0.218 


0.228 


.473 


.294 


.294 


.225 


.232 


.545 


.290 


.292 


.2.33 


.235 


.(U5 


.286 


.289 


.236 


.237 


.689 


.282 


.286 


.238 


.239 



HEAT AND COMPRESSIBILITY EFFECTS IN INTERNAL FLOW SYSTEMS 



33 



1.4 



1.2 



1.0 





























— 4 


\ 














J 




oer 


48. 








T. / 


f-r.. 




/ 


















Mo- 








































/ 




































+• — 
































































J 
















-CL. - 












_Xi 


*\ 










/ 






























\ 
\ 










t 








- /nrerna/ pressure coerr/u/en/ 
based on theoreiica/ analysis 


.A 








1 












o 


Ex 


per 


-/>77 


ent 


ol f. 


Doir 


7fS 






\ 








/ 


































\ 

\ 








































\ 




















Ve 






ra 


fio 




















































y 



















































.1 



.5 
l/L 



.6 



.8 



1.0 



Figure 18.— Tomperaturc ratio, internal pressure coefficient, and velocity ratio for tail B 
at A/=().(il5. Heat input. rKS.2 Btu per so.uare foot per second. 

This relation is shown for comparison with the oxp(M'ini(Mital 
data obtained from the wake pressure survc^ys in ligiire 10. 
Virtually all of the heat energy theoretically recoverable was 
actually recovered. The data shown re])r(^sent a very wide 
range of heating conditions and internal mass-flow rates. 
The maximum air-temperature rise in the tests was about 
:300'' F Mild tli(^ flow coefricicMit C vnri(Ml from about 0.025 
to 0.070. 

The recovery process was not measurably afl'ected by any 
increase in ducting loss chie to heating bcH'ause of the ex- 
tremely low initial value of the ducting loss for the mo(h4 
tested. With tail B, the medium-size outlet, the loss across 
the racHator was only 0.008(/o and the loss in the converging 
duct behind the radiator was about 0.007go for a mass-flow 
coefficient of 0.050. 

The heat loss through the walls between th<' radiator 
and the tail outlet is estunated as approximately 1 percent 
of the heat added. The heat loss determiiKMl IVoin tbe wall- 
temperature and bouiuhiry-layer surveys agreed with the 
theoretical calculations of this loss. Allowance for this 
small correction was madc^ in llic !'c(hi('ti()ii of (lie data, but 
there was no measurabh* cliaiigc in tlie rcMoveiy ])rocess with 
so small a heat loss. 

The force measurcincnls of th(^ iK^t thrust caused by 
heating pi-ovidc^l w means of (lel(M-inining a vahie of e^/ 
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that includes the external losses associated with exhausting 
the hot air of increased velocity. This value of €// is shown 
in figure 17. As noted from the definition of (see equa- 
tions (31) and (33)), the over-aU (^fliciency given by the 
force tests was reduced by the mechanical or Froude effi- 
ciency factor €m in computing €//. The force-test data 
closely approach the ideal value of e//. The results are con- 
sistent, however, in indicating small losses. Comparison of 
figures 16 and 17 leads to the conclusion that the losses 
must be due to slight increases in external drag from the 
heated exhaust flow. 

TABLE IV.— COMPARISON OF EXPERIMENTAL AND 
THEORKTICAL VALUES OF MASS-FLOW RATIO 
[Tail B, baffle plate] 



M 


Mass-flow ratio, pmlpa 




//»90 Btu/sec 


Equation 
(27) 


Equation 
(26) 


Experi- 
mental 


Equation 
(27) 


Equation 
(26) 


Exp(»ri- 
mcntal 


0. 263 
.347 
.403 
.473 
.545 
.615 
. (iSW 


0. 767 
.767 
.767 


0. 776 
.786 
.790 


0. 772 
.78! 
.786 














0. 621 
.6:39 
.657 
.668 
.676 


0. 675 
.69(i 

■ .719 
.734 
.751 


0. 672 
.694 
.718 
.736 
.756 








.767 


.828 


.819 



Measurements of the external (h'ag with the tail rake were 
made in only the vertical plane (flg. 14). The results so 
obtained showed no consistent external-drag changes with 
heating. In any case, such changes in external (h'ag w(M'(* 
very small; a 1-percent change in the ch-ag of the body would 
account for the small divergence between the ideal and the 
measured efficiency at the higlu^st test speed and heat input. 

The heated flow from the annular outlet had as little 
effect on the external drag as did that from the tail outlets. 
This result does not seem surprising when it is considercul 
that the dynamic pressure in the heated flow was virtually 
equal to that of the cohl outh't flow. The increased velocity 
is offset by the (l(M'i-(»as(Ml density so that thc^'e is little change 

in- P4V4^ and coircspondingly little change in external (h'ag. 
A 

(See i)art 1.) 

The efliciencies obtained with tlie chict pressure re(hiced 
])V 0.27go by means of the baflh' ])late approached the 
precHcted efticiency for the nieasunnl compression ratio as 
closely as in tlu* runs with low internal losses. 

Figures Hi and 17 indicate that the use of the ideal rela- 
tion for the heat-cycle efficiency is justified for design calcula- 
tions. If the duct-flow Alach number is low, say less 
than 0.10, the heat-cycle efficiency is adequately given by 
the simple relation of eciuation (35), 



Figure 19.— Variation of intomnl mass flow with outlot donsity ratio. 



Otherwise, the more exact relation (equation (o7)) nuist be 
'employed, in which case en may be obtained for the particular 

values of iXU and ^^^^ involved. 
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Tho data accumulat(Ml in tlic tests ar(» uscrnl in vcrifving 
the relations, dovolopcMl in part I, coucrrniug tho calculation 
of the flow characteristics at key stations in the chict, the 
effects of heat on the over-all flow characteristics, and the 
net propulsive force. Comparison of the computed and the 
measured values of these characteristics will now he made 
in the fore.u'oin^' ordiM*. 

Comparison of calculated and measured internal-flow 
characteristics. — ^The internal-llow eliaractcM'istics at the k(\v 
stations in the duct W(M-e ('()n)j)nte(l I'loni the equations of 
part I on the hasis of the lollowino; constants, which 
correspond to one of tlu* test conditions with tail B: 



//. I^t II per second. ()8.2 

r,„ 

miles jHT hour 456 

feet per second. 069 

po, pounds per square foot 1653 

Po, slug per cubic foot 0.001070 

Cor^ 75 

(A/?/(?o),/„^,, including dilfuser loss 0.03 

rjd 0.93 

I\ 0.11 

(outlet area), scjuare foot 0.0547 



The computed mass-flow coeflicirnt is 0.0399, as compared 
with the measured value, 0.0401. The area, .lo=0.0408 
square foot, which corresponds to the computed mass-flow 
coefficient, was used in calculating the j)ressure, the tempera- 
ture, and the velocity at each station in the duct according 
to the ecpiations of j)art I. The results are compared with 
the measured data in figure 18 and show that the flow 
characteristics of an effici(Mit syst(*m, in which tlie vcdocity 
distribution is fairly unifoiin. can he comput(Ml with good 
accuracy. 

Effects of heat and compressibility on the internal total- 
pressure loss. — Table III gives a comparison of tlu* nu^as- 
ured variations of total-|)r(»ssure loss for the tests with the 
internal resistance plate with the results comput(Ml by 
(Ujuation (23). Two rates of heating and a range of values of 
Mach numlxM- an* shown. The reduction of the t()tal-pr(»ssure 
loss with inci-eascul heating is du(^ to th(» (l(H*r(^ase in internal 
mass flow accompanying the addition of heal. For a given 
rate of heating, the total-])ressure loss dcH-reases as the 
Mach number incnnises in systems with a fixed outl(*t 
size. This result is caused by tlu^ a|)preciable inci-ease in 
denspity of the air at the radiatoi' section of the duct as the 
Mach number increases. AVith a given outlet opening and 
fixed rate of heating, the mass-flow ('(x^fficient doc^s not 
change greatly as the Mach numlxM* advances; {he increased 
density at the radiator section th(M"(^foi"(* rc^sulted in decreased 
velocity ratios through the radiator and in correspondingly 
decreased radiator losses. It wUl be notcul from table III 
that equation (23) predicts both the effects of heating and 
of Mach nunibei" veiw accurately. 

Effects of heat and compressibility on the internal mass- 
flow rate. — As shown by e([uation (2(V), tlu^ mass flow 
through a system with a given outlet ari^a will he i-educed 
in proportion to the squan* root of the outlet diMisity, 
provid(M[ that tlu^ intcM-nal losses are n(^gligibl(\ Figure 19 
comparers this th(M)retical variation witli the experinuMital 



results for all of the tests except the runs with tlu^ added 
internal resistance*. The agnn^ment is exc(»llent. For these 
tests, the decrease in mass-dow rate was independent of 
Mach number for a Mach number range from 0.2() to 0.75. 

WluMi an appreciable internal total-pressure loss is pn^sent, 
as in most practical installations, the change in mass-dow 
rat(* no longer varies simply as p^^^^ but depends on both the 
j internal resistance and on the Mach number Mq. Table 1\' 
shows the How variation obtaincul experimentally with the 
resistance plate iu the duct and, for comparison, the flow 
variations computed by equation (26) and also by tin* 
approximate e((uation (27). The exact equation (equation (26)) 
agrees very closely with the (*xperim(*ntal results, but 
equation (27) departs apprcn-iably from the experimental 
values at the higher Mach numbers. It has been found, 
however, that the approximate equation (27), which does 
not account for Mach nund)ei- chang(»s. yields sufficicMitly 
accurate results at Mach nund)ers below about 0.30. 

With an appreciable internal loss, the reduction in mass- 
flow rate with heating is less than the reduction indicated 
by the p^^^^ relation. The reduced total-pressun* loss accom- 
panying the addition of heat in a system with a fixed outlet 
opening tends to increase the internal flow and this ci\ort 
nullifies to some extent the throttling eftect of the lu^at. 

Computation of the net force due to the heated internal 
flow. An example of the simple method described in part J 
for computing the net force wdl now be given. The same 
duct characteristics and the same initial conditions assumcul 
in th(* preparation of figure 18 will be used. In addition to 
the (low characteristics previously listed, the total-pressure 

loss —=^0.238 and the coni|)ression ratio = 1 .22() are vc- 

(|uired. These valuers were computiMl from the basic data. 

Thv low-spe(Ml drag coefficient without lu^at is (from 
equation (39)) 

t7;;/ = 2X0.0401 (1-/1 -0.238) =0.0103 

The drag coeflicient without heat but at the design .Mach 
nundxM" .\/„ = 0.()15 is (from fig. 9) 

^7// = 0.0 103X0.790 = 0.0081 

From equation (35) (or fig. 8 for J/o~0), the efficiency of 
i the lu^at cycle is 

e//=l- (1.226) " -^•' = 0.0565 

Th(^ di'ag reduction due to heating is given by equa- 
tion (40). The value of e^f required in this formula is expres- 
sed by equation (34), which cannot be evaluated imtil 
equation (40) is cvaluat(Ml. As a trial value, e.u will be as- 
sumed (Hpial to 1.0. The drag reduction due to heating is 
tluMi, by equation (40) 

Kr JHe„eM _ 77S X 68.2 X 0.0565 X 1 .0 

<j,FVo ~ ' 441 X 1.009X669 - -^^-^^l^^l 

The \ alu(> of the net force is, finally, 

Cn= Cn + AC„ = 0.0081 - 0.0 101 = - 0.0020 
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11" this (ii-sl a|)j)i-().\iinatioii is used for the net foree, n vahie 
for e.u may he eompiit(^d from equation (M)] that is, 



1 



1- 



- 0.0020 + 0^0081 Y 
4X0.0401 / 



1.04 



With this vahie for e.^/, [he new value o))taine(l for the thrust 
(hie to h(»ating is 

ACn=-O.Om X 1.04= -0.0105 ' 

and the net foi'ce is 

C;> = 0.0081 -0.0105= -().()()L>4 

No further trials need he ma(h^ in this exam])le l)eeause e^f is 
so close to unity. The measured value of the net inteinal 
drag for this test condition was —0.0022. The slight dilfer- 
once between the measured and the computed results indi- 
cates the small error involved in the assumption of uniform 
outlet-velocity distribution. Th(^ foi'c^going UK^thod is exact 
except for this assumption. 

Comparison of subsonic ram-jet propulsion system with 
internal-combustion engine. The e()ni|)ressi()n ratio eorres- 
])on(ling to stagnation pressure^ in the duct at a Mach number 
of 1.0 is smaller than the values employed in internal- 
combustion engines. The ideal lu^at-cycle (^fHciency for the 
duct system considered in this report (ram j(4) is, therefore, 
lower at su])sonic speeds than for the internal-combustion 
(Migine. A number of energy losses, how(^v(M-, occui- in tlu^ 
use of the intiM-nal-eombustion engine that do not occur in 
th(* heat(H] duct (low. Thr^o losses may be lishul as follows: 

1. Mechanical losses in the moving j)ai1s of the (Migine 

2. Pow(M' i-equir(Ml foi- cooling-air flow 

3. Increase^ in aii'piani^ drag due to {\\r engine installa- 

tion 

4. Fi-iction and compressihility losses in the |)r()pell(M- 
Ki-om a considei-ation of these factors, a value of 24 ])ercent 
ioi- th(^ tluM-mal cHicicMicy of a typical (Migin(» installation w^as 
ai-i'ived at. This value does not include ihc induced losses 
in the propeller slipstream or the compressibility losses. In 
the cas(» of the duct flow, the present I'esults indicate that the 
ideal thermal efTiciency can closcdy approached if the 
ducting losses can be kept low. Furthermore, the results 
of i-eferences 4, 5, and 7 show tliat installation of the duct 
with suitable inlet and outlet openings can be accomplished 
without increasing the airplane drag. A comparison of th(^ 
two systems of propulsion is shown in figure 20. The ideal 
cfhciency of ihc duct system becomes of comparable magni- 
nide to the thermal efficiency of the engine^ installation at 
Mach numbers approaching 1.0. 

Tlu^ power recoverable from the cylinder-cooling How 
(^'Meredith effect," reference 1) in a typical air-cooled 
(Migine is also shown in figure 20. The pressure-drop results 
given in figure 6 were used to compute the compression 
ratio P2/P0' The heat dissipated in cylinder cooling was 
assumed to be 0.4 brake horsepower. As shown in the 
figure, at a given Mach number, greater recovery is possible 
at sea level than at altitude because the required Ap/qo in- 
creases rapidly witli increase in altitude. The energy 
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Kkji re 20.— romparison of the ovor-all oiricioncy when tho available heat is added iti an 
internal-combustion engine and in the duct at stagnation pressure. 

recovery for this air-cooled engine installation thus amounts 
to about 3 percent of the brake horsepower at sea l(^v(4 at 
A/o = 0.()0, but at 30,000 feet ther(^ would be no recovery at 
this speed. For Iow(m- si)ee(ls at 30,000 feet the efiect of the 
heat would cause a j)ow(M- loss. 

SUMMARY OF RESULTS 

COMPRESSIBILITY EFFECTS 

The n^sults of both the analytical and experiuKMital ])arts 
of this investigation are summarized as follows: 

1. Compression of the air in th(> duct ahead of the radiator 
results in largc^ increase's in pressure, temperature, and 
density at the flight Mach numbers now connnonly attained. 

2. The Mach number of the flow in the radiator is an 
important parameter governing the flow changes across 
the radiator or air-cooled engine. At high internal Mach 
numbers, the density decrease through the radiator due to 
the heat added is greatly augmented by the effects of friction 
and momentum change. 

3. A large increase in the pressure drop, corresponding to 
this density decrease, occurs for a given mass flow^ as the 
Mach number advances. This increased pressure drop 
results partly from increased friction losses and partly from 
the change in momentum of the air flow within the radiator. 
The magnitudes of these compressibility or density-change 
effects increase with radiator drag coefficient as well as 
with the Mach number. 

4. In order to maintain required cooling of air-cooled 
engines as the altitude is increased, an increasing volume 
rate of air flow is required because of the rednction in atmos- 
pheric density. The Mach number of the internal flow 
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tlirroforo incroasos with altitudo and the compressibility 
effects consequently become more severe. It appears that 
the use of a blower to produce adequate pressure drops at 
high altitudes will b(* i-e(juired for some existing (engines. 

5. The maximum possible mass-flow rate is attained when 
sonic velocity is reached at any station in the internal flow 
system. Increasing the pressure drop after this limiting 
condition is reached will not increase the mass-flow rate; 
the added energy will be dissipated in internal shock waves. 

6. The attainment of sonic velocity is the limiting possible 
flow condition regardless of the rate of heating or the mag- 
nitude of the internal-friction losses. Increases in either 
of these factors, how(*v(M-, reduce tlu^ mass-flow rat(* obtained 
at the limiting flow condition. 

7. In conventional cooling-systcMu ariangemeiits, sonic 
velocity would first be obtained either within the radiator 
at the tube exits or at the outlet opening. Existing air- 
cooled engines generally employ much higher cooling-air 
velocities than an efficient radiator, and the critical flow 
condition might occur under flight conditions now attainable. 
Calculations for a typical air-cooled-engine installation 
indicated that sonic velocity would be attaincMl in th(» baffle 
exits at altitudes ranging from about 40, ()()() to 45,000 feet 
if a maximum head temperature of 450° F and a fuel-air 
ratio of 0.10 were maintained. 

EFFECTS OF HEATING 

1. In usual internal flow systems heat is added in an 
expanded section of the duct where the static ^Dressure is greater 
than atmospheric pressure. Under these conditions the 
momentum and the kinetic energy of the flow in the wake 
are increased. The nit'io of the increased kinetic energy to 
the heat energy addcnl is denoted as the efficiency of the 
energy-recovery process. The momentum increase, which 
determines the thrust or drag reduction due to heating, is 
related to the kinetic-energy gain by a factor depending on 
the relative magnitudes of tlu^ wak(^ velocity and the flight 
speed. 

2. The tluM'mal efficiency of the recovery process depends 
mainly on the ratio of the static pressure in the duct at the 
entrance to the heating chamber to the atmospheric pressure 
and, to a lesser extent, on the Mach number of the flow in 
the duct, and on the rate of heating. Increases in either the 
Mach number of the duct flow or the heating rate cause 
decreases in the efficiency. 

3. An important effect of addition of heat in a constant- 
area duct is to increase the Mach number of the flow. For 
heating rates corresponding to temperature rises of the order 
of 2500° F, a Mach luunbcM' of 1.0 is reached at the end of 
the duct for an inlet Mach number of about 0.2. This criti- 
cal or choking inlet Mach number dictates the area expan- 
sion required in th(* diffuser ahead of the combustion chamber. 

4. The compn^ssion ratio, and hence the efficiency, in- 
creases approximately as the square of the flight Mach num- 
ber. Maximum recovery at a given flight Mach number 
requires that factors tending to reduce the compression ratio 
be minimized. Poor installations having large internal- 
pressure drops will suffer also from decreased gain from the 
heat-recovery process. 



5. The power recoverable from the heat added to the cool- 
ing air (Meredith effect) was found to be about 3 percent of 
the brake horsepower for an existing air-cooh^l engine at a 
flight Mach number of 0.60 at sea level. At 30,000 feet, 
however, no recovery occurred at this flight Mach number, 
and at lower Mach numbers the recovery was negative. 

6. The addition of heat in an internal flow system lias a 
throttling eft'ect on the mass-flow rate. The flow for a given 
outlet area decreases approximately in proportion to the 
square root of the outlet density. It is essential that this 
effect be allowed for in computing the required outlet- 
opening size. 

7. The outlet dy nam ie pressui'e ~ p^Vi' does not change 

appreciably with heating because, in most cases, the reduc- 
tion in density is nearly offset by an increase in tlu^ outlet 
velocity term ^4". The drag of bodies in the wake, therefore, 
would not be changed due to heating, except for possible 
scale effects associated with tlu* r(»duc(Hl Reynolds number 
in the heated flow. 

CHARTS AND EQUATIONS 

For eonvcMiience, the formulas and charts showing the 
princii)al lieating and eompr(»ssibility effects an* listen! as 
follows: 

(a) Pressure, dcMisity, and temperatui-e ahead of radiator: 
figures 2, 3, and 4, and equation (5) 

(b) Density decrease within the radiator and across the 
radiator: figure 5 

(c) Pressure drop within the radiatoi*: (Hjuation (11) with 
the aid of figure 5 

(d) Pressure drop across radiator: efjuation (13) or (14) 
with th(* aid of figure 2 and figure 5 

((») Effect of d(Misity change across radiator on the pres- 
sure di'op: cfjuation (14) 

(f) T(Mnperatur(» rise across radiator: equation (18) 

(g) Required outlet area: equation (28) with the mass- 
flow ratio given by equation (26) or (27) or by figure 7 

(h) Ideal efficiency of heat-recovery cycle: equation (35) 
or (37) or figure 8 

(i) Thrust due to the addition of lunit: e(|uation (40) 

(j) Net force due to luxated inteinal flow: equation (41) 
with the aid of figure 8 

CONCLUSIONS 

1. Th(» flow details in an efficient internal flow system with 
addition of heat can be computed with good accuracy by the 
one-dimensional compressible flow relations presented in 
part I. 

2. The actual thermal efficiency obtainable from a ram-jet 
type of propulsion system is virtually equal to the ideal 
thermal efficiency. The thrust available from ram-jets can 
thus be calculated with good accuracy from simple theory 
if the rate of heat addition to the flow is known accurately. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., July 21, 1942. 



APPENDIX A 



CALCULATION OF PRESSURE DROPS REQUIRED FOR COOLING A TYPICAL AIR-COOLED ENGINE AT HIGH ALTITUDES 



Basic experimental data. — Test data used to establish the 
basic sea-level air requirement and pressure drop were 
obtained on a Pratt & Whitney R-2800-27 engine installed 
in a nacelle in the Langloy 16-foot high-speed wind tunnel. 
The operating condition mikI the basic data obtaiiKul were as 
follows: 

Brake horsepower at 2400 rpiii (normal rating), assinned con- 
stant with altitude 1600 

Fuel-air ratio 0.10 

Air temperature, T2 

OF 100 

°F abs . 560 

Maximum head temperature 

op 450 

°Fab8 910 

Average head temperature 

OF 405 

°F abs 865 

Mass flow of cooling air 

slugs/sec 1.165 

lb/sec 37.5 

Pressure drop Apr, in. of water 15 

Frontal area A2, sq ft 17 

Kquivalent leak area ratio ^r/^U' ^-14 

Variation of mass flow of cooling air with altitude.— ^The 
mass [low recjiiiicd to maintain a fixed head temperature 
depends on a number of variables, but a discussion of the 
details of the variation with altitude is considered beyond tlie 
scope of this paper. A relation currently in use for deter- 
mining the mass flow will therefore be used without explana- 
tion of its derivation. If a supercharged engine with inter- 
cooler is assumed, this relation for the mass-flow variation is 

X( ^ head sea leve l 

.„ (T^ T \^-^ 

head •*■ 2) altitude 

mass-flow calculations will be cai-i i(Ml out for values of 
the fliglit Macli number Mo ecjual to 0.4, 0.6, and 0.8 through 
tlu^ altitude range. As pointed out in part I, the tempera- 
ture abead of tlu^ (Migine To depends on the flight speed as 
well as on the atmos])lieric t(Mn])erature. For Army Summer 
Air, the value of T2 was detcM mined from figure 4. For Mq 
(Mjual to 0.6, the mass flows computed at altitude ])y the 
foregoing relation ai*e as follows: 



Altitude 


To 




Mass flow 


(ft) 


(°F abs.) 


rr'abs.) 


(lb/sec) 


0 


560 


000 


46.3 


10,000 


524 


5t)2 


38. 0 


20.000 


489 


524 


31.7 


30,000 


453 


480 


27.0 


40,000 


418 


448 


23.4 


50, 000 


400 


429 


21.9 



The im])ortance of allowing for tlie tempei-alure rise caused 
by adiabatic compression ahead of the engine is illustrated 
by the sea-level case for which the mass flow was only 87.5 
pounds per second, the compression eftVct being neglected, 
as com])ar(Ml witb 46.3 ])oun(ls ])er scm'OihI shown in the table. 

Details of the pressure-drop calculations. — The required 
])ressure drop will be computed as the dro]) between stations 
2 and v., plus the drop in the bafflers. If (ujuation (13) is to 
be applied to an air-cooled cylinder on which the baffling 
covers only the rear half of the cylinder, account must be 
taken of the fact that approximately half of the heat is 
added to the cooling air ahead of the baflle at virtually 
constant pressure. Thus, the density of the air immediately 
ahead of the baffle at station 2a is obtain(Ml from tlu* reflation 



P2a 



H/2 



CpgmT2 



The value of H, the IhmU added to th(^ cooling air, was taken 
as the heat (^quivalcMit of 0.4 brake* hors(*power. Thus, 
equation (13) for the air-cooled engines becomes 



(l2 



-i^K-i:^*)+p7:(t)'K-<?::-0] 



in which Prjp-i and poalPr,, are obtained from figures 2 and 3, 
respectively, for the design valuers of Mr^ and ArjAo. 

The ratio Prjpr^ is obtainc^l directly from figure 5 for the 
design values of M,^ and Co^. The friction component of the 
engine drag Cd^ is assumed to result only from friction within 
the baffles and was obtained from tlu* basic engine data as 
follows: 

(1) The over-aU test pressure drop Apr was 78 Ib/sq ft 
and P2a/P3 was 1.08. 

(2) Substituting tlu^se values in equation (14) gives 

7S = Apf,X^ (1 + 1.08)4-2X1.07X51 (1.08-1) 
+ 1.07 (51-1) 



or 



(3) 



Ap,=l5.2 Ib/sq ft 
^P/i 15.2 



''friction 

(4) By equation (8), 

_ 14.2 22 
^^/~ 2 



1.07 



14.2 



^l + ^'O 

r\ PrJ 
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As ail oxamplo, tho results of tho calciilatious for a typical 
iVi^ht speed and altitude are listed in the following table: 

[Altitude 35,000 ft; A/^d=0.60] 



Station 
(See 
fig. 1.) 


V 

Vsq ft/ 


p 

/slug \ 
Vcu ft / 


T 
(O p) 




M 


q 

Vsq ft/ 


O 


497 


0. 000666 


435 


0. 60 


125 


2 


630 


. 000788 


466 


.0545 


1.31 


20 


630 


. 000730 


503 


.0567 




r2 


548 


. 000637 


483 


.456 


79.5 


rz 


498 


. 000543 


514 




.518 





The pressure (h'O]) across the engine^ is 

A2;r=2^2— i>r3=l32 Ib/sq ft 

See figure 6 for pressure drops for other altitudes. 

In order to illustrate the large error involved in the 
neglect of the effects of the density decrease across the engin(\ 
the pressure drop will be computed without regard for these 
eflfects. In this case, the usual assumption is 

— '"=Constant with altitude 

<l2 



or 

- XA2?rX V-=Constant with altitude 

wluMice 

^^cn leiel\ nrfl . . / n . 

\''''sea lerel/ Paltilutle 

= 94 Ib/sq ft 

as compared with the exact value, 182 Ib/sq ft. (See fig. 6(b).) 

Th(* api)roximate method (equation (14)) for obtaining 
the pressure drop across the engine without detailed calcu- 
lation of the flow within the bafHc^s was found to yield 
approximately the same* results as the foregoing more exact 
analysis except when the velocity of sound was approached 
in the baffles. The approximate method gives no indication 
of the nttninnient of sonic velocity. 



APPENDIX B 



REDUCTION OF DATA FOR TESTS WITH HEAT ADDITION 
INTERNAL DRAG 



No information on tho snbjoct of compiitini^-, from wako- 
survey data, tho drag due to hoatod iiit(Miial flow at high 
speeds has been pubhshed. Th(* nu^lhod to be described is 
op])licabl(' to the usual oxi)(M*imental case whore the rate of 
healing and tlie static and total pressures at the outlet are 
measured. The momentum relation (equation (21)) is the 
basis of tho calculations. ^Mien allowance is made for the 
nonuniformity of the flow, equation (21) becomes 



or 



It is convoniont to wi-it(^ this rolation in tlu^ [ovm 



dA (42) 



If the usual assumption is made that the total pressure at 
station 5 (wIhm'o 2>5=2^o) is the same as the outlet total 
pressure, tho vai ious terms in this relation may be evaluated 
as follows: From tho isonli'o])io relation, 



and, also, 



1/2 



\k-vJ L(i+^)J 



in which llio total and static pressures are measured and 
l + rj is tho compressibility factor. A more useful formula 
for 1+t; than that given in tho list of symbols involves tho 
total and slat ic ])rossures as measured in the wake ratluM' than 
tlio wak(^ Macli nunilxM'. This r(^Iation is as follows: 



l+,= l+().357 ^;;^^;^ + «.05,[^f,4yJ 



1+ 



Similarly, since h. 

\qo/ \ho-poJ L(1+^)5J 

The only remaining t(Min r(Hiuirod for the evaluation of 
oquation (42) is (po/po)'^'. Tho wak(^ density ratio pjp,, 



depends on two factors: (1) the rate at which heat is adch^l 
to the internal flow; and (2) tho amount of aerodynamic 
heating, which in turn is govorn(*d by the Mach number 
AIq and by tlu^ internal loss. Tlw term (ps/po)^^^ may bo 
evaluated by writing tho energy equation between stations 0 
and 5. 

— +-77- +^^7^*.-/o+ -^— + -:y+JgCvl5 
Po m ps z 



R 



and Mo^ 



whonc(v 



jc ivlations p5=7>o, TQ=T^{prJ p^) , po = Po(lRTo, c,= j^^_^y 
PqVq^ 



are substitut(Ml in this (Kpiation, then 



P5 



2^ Po 

(lo pi 



J^CpgtnTo 



1/2 



l + 0.20Mo\qM 



l+().20Mo'+ 



Cp(jmTo__ 



1/2 



Thoso relations for {pjp.y^', (q^qoY^', {qM'". and {pjp,yi\ 
when substituted into equation (42), determine the internal 
drag corresponding to tho point in tlu* wake at which tho 
pressures are measured. Integration across the outlet open- 
ing gives the required value of the total drag coefficient due 
to tho int(n-nal flow. 

EXTERNAL I)RA(; 

The equations for tho internal (h-ag apj)ly for tho external 
drag with tlu^ follow ing changes: 



(1) 



P,y ri +0.20Mo^(g5/go) "| 
~L l+0.20Mo2 J 



(2) The integral focpiation (42)) is ovahiatod botw(HMX the 
surface of tho mo(U'l arid tlu^ outer extrennly of the wake. 

INTERNAL MASS-FLOW COEFFICIENT 



C- 



PqFV 



(43) 



These are the same quantities as already derived in connection 
with the internal drag cooflicient. 

CALCULATION OF HEAT-CYCLE EFFICIENCY 

The heat-cycle efficiency (equation(31)) can be very con- 
veniently computed from the difference between the drag 
ooeflPicient for the cold condition and the drag coefficient for 
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th(' boated condition. From equations (33) and (34), 
" eM JH^M J II \ e^f J 

=«gl"j(0.'-c„)[i-4<^]] ,«) 

Inasmueh as the mass flow is less for ihc hnvUnl coiidilion 
than for the cold condition, the measured internal (h-ag 
without heat is therefore greater than tlie dc^sinnl valu(» D' 
which would be measured if the mass flow (hd not change. 
For smafl variations in mass flow, however, it is possible to 
correct accurately the (h'ag for the measured unheat(Ml coiuh- 
tion to the hypothetical unheated condition. The relation 
between the total-pressure loss for the desired unheated con- 
dition and for the measured unheated condition is 

The outlet total pressure foi- the desired condition is, then, 




With this value for /14, the desir(Ml value of Co' is computed 



directly from ecpiation (42). This value of Ci/ permits 
computation of from ec[uation (44). 

Equation (44) may be used to evaluate either for the 
wake-survey values of Cd and Co or for the force-test values 
of thes(» drag coefficients. 
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I 




Y 

Z 

Positive directions of axes and angles (forces and moments) are shown by arrows 



Axis 


Force 
(parallel 
to axis) 
symbol 


Moment about axis 


Angle 


Velocities 


Designation 


Sym- 
bol 


Designation 


Sym- 
bol 


Positive 
direction 


Designa- 
tion 


Sym- 
bol 


Linear 
(compo- 
nent along 
axis) 


Angular 


Longitudinal 


X 


X 


Rolling 


L 


y — >z 


Roll 




u 


P 


Lateral- ___ 


Y 


Y 


Pitching 


M 


z — >x 


Pitch. 


e 


V 


Q 


Normal 


Z 


Z 


Yawing 


N 


X — >y 


Yaw - 




w 


r 





















Angle of set of control surface (relative to neutral 
position), 8. (Indicate surface by proper subscript.) 



4. PROPELLER SYMBOLS 



D 

V 


Diameter 
Geometric pitch 




P 


P 

Power, absolute coefficient Cp= — 


p/D 
V 


Pitch ratio 
Inflow velocity 






Speed-power coefficient= ^ 


V. 


Slipstream velocity 




V 


Efficiency 


T 

Q 


Thrust, absolute coefficient Ct= 
Torque, absolute coeflacient Cq- 


T 
Q 


n 


Revolutions per second, rps 
Effective helix angle— tan ^^27rrn) 



5. NUMERICAL RELATIONS 

1 hp=76.04 kg-m/s-550 ft-lb/sec 1 lb=0.4536 kg 

1 metric horsepower =0.9863 hp 1 kg=2.2046 lb 

1 mph=0.4470 mps 1 mi= 1,609.35 m=5,280 ft 

1 mps=2.2369 mph 1 m= 3.2808 ft 



Absolute coefficients of moment 

^'~qbS ^''~qcS ^""gfiS 

(rolling) (pitching) (yawing) 



